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Abstract 

During  June  1988,  an  experiment  using  NORDA's  towed 
underwater  pumping  system  (TUPS)  was  conducted  in  shallow 
water  off  Panama  City,  Florida  to  gather  data  in  support  of 
the  Navy's  Airborne  Bathymetric  Survey  System  (ABS)  . 
Continuous  measurements  were  made  of  upwelling  irradiance, 
fluorescence,  t ransmissometry,  water  depth,  temperature,  and 
salinity  using  sensors  in  TUPS.  Ambient  light  measurements 
and  navigation  data  were  recorded  from  instruments  aboard  the 
research  vessel.  This  report  describes  the  experiment,  the 
cruise  itinerary,  and  the  methods  used  to  carry  out  out  the 
measurements.  Data  are  presented  in  graphical  form,  along 
with  a  description  of  the  data  reduction  protocols  and  some 
preliminary  results.  Listings  of  computer  programs  used  for 
processing  the  data  are  included.  A  description  is  given  of 
the  contents  of  a  magnetic  tape  (provided  to  NORDA) 
containing  all  programs,  raw  and  processed  data. 

Preliminary  analysis  of  the  data  indicate  that  the 
quality  of  the  light  sensor  and  depth  data  are  good.  A 
significant  inverse  relationship  between  depth  and  upwelled 
irradiance  was  observed  under  some  conditions.  The 
relationship  holds  well  in  areas  of  bright  sand,  but  is  more 
complicated  in  areas  of  variable  bottom  type  and  water 
clarity.  The  broad  suite  of  light  and  water  quality  sensors 
aboard  TUPS  make  it  an  ideal  system  for  providing  the  data  to 
resolve  such  problems . 
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1.0  INTRODUCTION 


1.1  Rationale  for  the  Study 

The  Naval  Ocean  Research  and  Development  Activity 
(NORDA)  has  been  developing  and  testing  an  Airborne 
Bathymetric  Survey  System  (ABS)  to  measure  water  depth  in 
coastal  areas.  The  ABS  is  being  developed  by  NORDA  for  the 
1  Defense  Mapping  Agency  (DMA)  and  the  Oceanographer  of  the 

Navy  (CNO  OP-096)  under  the  NORDA  Coastal  Hydrographic 
Techniques  Program.  The  ABS  uses  a  laser  system  in 
•  conjunction  with  a  multispectral  scanner  which  measures 

upwelled  irradiance  at  several  wavelengths.  The  laser  depth 
measurements  and  the  multispectral  scanner  data  are  merged 

►  together  to  provide  an  algorithm  which  allows  computation  of 
depth  from  the  multi-spectral  scanner  data  alone. 

The  goal  of  the  work'  reported  here  was  to  gather 

>  environmental  data  for  making  ground-truth  measurement  of  the 
ABS.  In  addition  to  the  ground  truth  support,  this  work  was 
undertaken  to  help  improve  the  present  algorithms  used  by  the 

)  ABS,  by  gaining  a  better  understanding  of  the  ocean's 

optical-environmental  parameters.  To  accomplish  these  goals, 
NORDA  equipped  its  towed  underwater  pumping  system  (TUPS) 

I  with  a  depth  sensor  (0-100  meter  capable  with  2  cm 

resolution)  and  upwelling  irradiance  sensors  with  fixed 
wavelength  filters  of  465  nm,  507  nm  and  532  nm.  This 
)  equipment  was  installed  in  TUPS  during  the  spring  of  1988. 
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It  was  decided  to  perform  a  field  experiment  using  TUPS  to 
provide  ground-truth  coverage  for  ABS  overflights  to  be 
conducted  during  June  1988  off  Panama  City,  Florida  in  the 
Gulf  of  Mexico. 

For  this  experiment,  the  Geochemical  and  Environmental 
Research  Group  (GERG)  of  the  Department  of  Oceanography  at 
Texas  A&M  University  was  contracted  (ONR  Contract  No.  N00014- 
88-K-6003)  to  provide  the  ship  and  personnel  required  to  tow 
the  NORDA  TUPS  off  Panama  City,  Florida.  This  effort 
included  chartering  a  suitable  vessel,  installing  the  TUPS 
hardware  and  computers  aboard  the  vessel,  operating  the 
system  at  sea  during  the  cruise,  collecting  the  TUPS  data 
along  with  navigation  and  appropriate  atmospheric  information 
while  at  sea,  and  processing  the  data  upon  return  to  the 
laboratory . 

This  report  describes  the  cruise  data  collected  with 
TUPS  off  Panama  City,  Florida,  on  June  19,  20,  21,  1988.  The 
data  analysis  procedures  are  described  in  detail  and  copies 
of  FORTRAN  computer  programs  use  to'  process  the  data  from 
collection  to  presentation  are  provided.  A  complete  set  of 
all  data  and  programs  is  available  on  a  VAX-compatible  9- 
track  magnetic  tape  which  was  supplied  to  the  NORDA 
Scientific  Officer,  Mr.  Robert  A.  Arnone,  along  with  this 
report.  The  complete  data  will  be  processed  further  by  NORDA 
scientists  to  better  understand  the  relationships  between 
light,  water  clarity  and  depth  in  coastal  regions. 


1.2  Reasons  for  Using  TUPS 


To  provide  the  ground-truth  data  for  the  Airborne 
Bathymetric  Survey  System  (ABS),  NORDA  chose  to  use  its  towed 
underwater  pumping  system  (TUPS) .  TUPS  was  designed  and 
constructed  at  NORDA  for  a  project  to  study  "Chemical 
Dynamics  in  Ocean  Frontal  Areas".  To  meet  the  needs  of  that 
project  and  to  provide  the  capability  of  performing  other 
surface-water  surveys,  the  TUPS  tow  vehicle  was  designed 
(Rein  et  al.,  1985)  to  carry  a  large  suite  of  oceanographic 
instruments.  It  was  also  designed  to  be  easy  to  modify  in 
order  to  allow  other  instruments  to  be  quickly  installed  and 
interfaced  to  its  on-board  computer.  The  instrument  layout 
of  TUPS,  as  originally  configured,  is  shown  in  Figure  1.  The 
tow  vehicle  configured  in  this  manner  had  been  used 
successfully  to  study  variability  in  water  optical  properties 
and  environmental  parameters  in  the  western  Mediterranean  Sea 
(Arnone  and  Wiesenburg,  1988)  .  Its  successful  use  in 
previous  studies  and  its  easy  adaptability  made  it  an  ideal 
instrument  for  studying  the  optical-depth  relationships  in 
coastal  waters. 


2.0  DATA  COLLECTION 


2 . 1  Study  Area 

The  area  chosen  by  NORDA  for  this  study  was 
area  south  of  Panama  City,  Florida.  This  site 
because  it  was  an  region  where  different  bottom 
available  in  close  proximity  to  each  other,  thus 


the  coastal 
was  chosen 
types  were 
minimizina 
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the  ship  travel  time  required  to  study  several  different 
areas.  i'he  Panama  City  site  was  also  selected  for  its  close 
proximity  to  KORDA  and  the  ability  of  the  ABS-equipped  P-3 
aircraft  to  operate  effectively  out  of  nearby  Eglin  Air  Force 
Base  . 

The  study  area  is  shown  in  Figure  2.  This  map  is  a  55% 
reduction  (original  scale  1:25,000)  of  a  section  of  KOAA 
National  Ocean  Survey  Chart  No.  11391  (DMA  Stock  No. 
11BHA11391) .  The  map  encompasses  the  area  where  the  TUPS  was 
towed  during  this  experiment.  On  June  19  and  20,  1988,  TUPS 
was  towed  in  the  Gulf  of  Mexico  south  of  Shell  Island.  On 
June  21,  1988  TUPS  was  .owed  in  the  shallow  area  north  of 
Shell  Island. 

2.2  Cruise  Itinerary 

The  cruises  undertaken  during  this  experiment  were 
conducted  aboard  the  vessel  Captain  Graydon  York,  a  forty- 
five  (45)  foot  crew  boat.  The  vessel  was  loaded  with  the 
TUPS  and  TESS  equipment  while  docked  at  Sun  Harbor  Marina  on 
June  18,  1988.  TUPS  was  configured  to  be  towed  from  the  port 
'ide  of  the  vessel  using  a  davit  that  had  been  constructed  at 
Texas  A&M  University  especially  for  this  experiment. 

The  Captain  Graydon  York  departed  Sun  Harbor  Marina  at 
1436Z  (0936  CDT)  on  June  19,  1988  and  proceeded  to  a  known 
reference  position  (30°0C.72'N,  85.41.60‘W)  to  calibrate  the 
TESS  LORAN-C  system.  The  LORAN  was  correct  to  0.01'  of 
latitude  and  longitude.  The  vessel  then  proceeded  out  the 


channel  and  south  of  Shell  Island. 


TUPS  was  launched  at 


1600Z  and  after  proceeding  seaward  was  towed  along  lines 
perpendicular  to  the  beach  to  provide  a  complete  survey  of 
the  area.  As  the  vessel  approached  the  beach,  we  endeavored 
to  go  into  as  shallow  water  as  possible.  On  two  occasions 
( 1 8 4 7 z  and  1854Z)  the  bow  of  the  boat  hit  the  bottom  and 
severe  course  adjustments  had  to  be  made. 

The  TUPS  survey  south  of  Shell  Island  was  completed  at 
2016Z.  The  TUPS  tow  vehicle  remained  in  the  water  and 
collected  data  while  anchor  Stations  1  ( 2 0 17 Z-2 1 02 Z , 
3  0°05 . 0  6N,  8  5°4  0 . 7  4 W)  and  anchor  Station  2  (2114Z-2202Z, 
30°04.56N,  85°41.05W)  were  occupied.  The  TUPS  tow  vehicle 
was  hoisted  back  aboard  the  vessel  at  the  termination  of 
anchor  Station  2  and  the  vessel  returned  to  port  at  2252Z. 

On  June  20,  1988,  the  Captain  Graydon  York  departed  its 
dock  at  1246Z  and  again  performed  a  LORAN-C  check  at  1305Z. 
LORAN  positions  were  exact  with  the  reference  position.  The 
early  portions  of  this  day  were  devoted  to  a  vertical  station 
with  the  NORDA  scanning  radiometer  system.  Station  3  was 
occupied  south  of  Shell  Island  at  30°02.15N,  85°45.95W  from 
1349Z  to  1440Z . 

After  leaving  this  station,  the  vessel  proceeded 
shoreward  to  launch  TUPS  in  shallow  water.  TUPS  was  launched 
at  1504Z  (1004  CDT)  and  a  zig-zag  pattern  was  run  toward 
Shell  Island  with  turns  at  1542Z,  1552Z,  1558Z,  1605Z,  1611Z, 
and  1615Z.  During  the  turn  at  1615Z,  the  power  generator 
aboard  the  Captain  Graydon  York  stopped  abruptly.  The 


generator  had  overheated  due  to  a  failure  of  the  water  pump 
on  its  cooling  system.  TUPS  was  retrieved  at  1645Z  and  the 
vessel  returned  to  its  dock  to  affect  repair  of  the  faulty 
water  pump.  No  further  data  were  collected  on  June  20,  1988. 
Also,  the  data  that  had  been  collected  were  unreachable  at 
that  point  as  neither  the  TUPS  or  TESS  data  files  had  been 
closed  properly  when  power  was  lost  to  the  computer.  The 
generator  was  repaired  during  the  evening  of  June  20,  1988. 

On  June  21,  1988,  the  vessel  left  port  at  1105Z  (0605 
CDT )  proceeded  to  the  LORAN  reference  point  and  then  to  a 
point  on  the  leeward  side  of  Shell  Tsland  (30°05.96N, 
85°41.31W)  where  TUPS  was  launched  at  1155Z.  TUPS  was  towed 
until  1430Z  when  the  vessel  anchored  for  Station  4  at 
30°06.27N,  85°41.72W.  Station  4  was  occupied  from  1430Z  to 
1507Z.  Station  5  (30°06.06N,  85°41.48W)  was  occupied  from 
1518Z  to  1552Z  and  Station  6  (30°05.85N,  85°41.36W)  was 
occupied  from  1602Z  to  1623Z.  At  this  station,  the  TUPS  tow 
vehicle  was  lowered  from  its  normal  tow  depth  beginning  at 
1602Z  down  to  near  the  bottom  and  then  returned  to  its  normal 
two  depth  at  1613Z.  Data  at  Station  6  were  collected  with 
TUPS  until  1621Z . 

The  Captain  Graydon  York  then  proceeded  to  Station  7 
(30°05.97N,  85°41.08’W)  which  it  occupied  from  1630Z  to 
1762Z.  At  this  station,  another  vertical  cast  was  made  with 
TUPS.  The  descent  started  from  the  surface  at  1643Z  and 
terminated  1.17  m  above  the  bottom  at  1651Z.  The  TUPS  tow 
body  was  retrieved  slowly  and  arrived  back  at  its  normal  1.0 


meter  tow  depth  at  1722Z.  After  the  end  of  Station  7,  the 
vessel  proceeded  to  Station  8  at  30°05.57'N,  85°40.61'W  which 
it  occupied  from  1740Z  to  1812Z.  At  the  end  of  this  station 
TUPS  was  retrieved  and  placed  in  its  cradle  in  order  to 
proceed  rapidly  to  a  point  outside  the  jetties  where  it  would 
be  launched  for  a  final  survey  through  the  Panama  City 
channel . 

TUPS  was  launched  again  at  1841Z  at  30°06.79'N, 
85°44.13'W  and  the  vessel  proceeded  on  a  course  of  320°T  at 
5.0  knots.  This  run  across  the  channel  was  conducted  from 
1841Z  to  1845Z.  The  transect  was  run  about  halfway  between 
channel  markers  1-2  and  3-4.  After  the  channel  crossing  the 
vessel  headed  south  then  into  the  channel  (at  1850Z)  south  of 
channel  markers  1-2 .  The  TUPS  was  then  towed  up  the  Panama 
City  channel  to  the  LORAN-C  reference  point  (30°08.72'N, 
85°41.60'W)  where  TUPS  was  retrieved  at  1926Z.  The  Captain 
Graydon  York  then  returned  to  its  dock,  arriving  at  1943Z 
(1443  CDT) .  The  NORDA  equipment  was  off-loaded  at  that  time 
ending  the  Panama  City  field  experiment. 

2.3  TUPS  Sensor  Configuration 

During  the  Panama  City  experiment,  the  towed  underwater 
pumping  system  (TUPS)  was  configured  slightly  differently 
than  shown  in  Figure  1.  The  original  configuration  included 
pitch  and  roll  sensors  (Trim  Pots),  had  only  the  capability 
for  two  upwelling  light  sensors  and  had  no  depth  measuring 
capability.  For  this  experiment,  the  trim  pots  were  removed 


to  provide  two  extra  data  channels.  One  was  used  for  a  light 
sensor  and  the  other  was  used  for  an  echo  sounder. 

A  list  of  sensors  used  on  the  TUPS  during  this 
experiment  is  shown  in  Table  1.  Three  upwelling  light 
sensors  at  fixed  wavelengths  {465  nm,  507  nm  and  532  nm)  were 
mounted  in  the  lower  quadrant  of  the  tow  vehicle  along  with 
an  Ulvertech,  Ltd.  Model  205  Echo  Sounder  (Depth  Sensor) . 
The  echo  sounder  produced  500  kHz  pulses  with  a  power  of  150 
watts  and  beam  width  of  four  degrees.  The  unit  produces  an 
output  voltage  of  0-10  VDC  which  is  proportional  to  0-100 
meters.  It  has  a  resolution  of  2  cm. 

Although  the  Ulvertech  echo  sounder  was  not  provided 
with  calibration  data,  calibrations  were  made  by  actual  depth 
measurements  (using  divers)  at  several  stations  during  the 
cruises.  The  echo  sounder  reported  within  1  cm  of  the 
measured  depth  at  3.1  meters  and  within  10  cm  of  the  reported 
depth  at  4.4  meters. 

2.4  TESS  Sensor  Configuration 

To  collect  incident  light  measurements  for  comparison 
with  the  TUPS  data,  we  used  NORDA's  TUPS  Environmental  Sensor 
System  (TESS)  .  This  system  (not  to  be  confused  with  the 
Tactical  Environmental  Satellite  System)  was  used  to  record 
total  irradiance  using  two  Eppley  pyroheliometers  as  sensors. 
These  measurements  are  important  since  clouds  passing 
overhead  reduce  the  amount  of  light  hitting  the  ocean  surface 
and  consequently  reduce  the  amount  of  upwelling  irradiance. 


Table  1.  TUPS  Sensor  Configuration,  June  1988 


Channel*  Sensor 


S/N  Calibration-*- 


FO 

Temperature  sensor. 

Sea  Bird,  Inc. 

SBE3-638 

9-18-87 

FI 

Conductivity  sensor, 

Sea  Bird,  Inc. 

SBE4-234 

9-18-87 

AO 

465  nm  Upwelling  Light, 
Biospherical  Instr. 

MCP-200H-7129 

5-2-88 

A1 

507  nm  Upwelling  Light, 
Biospherical  Instr. 

MCP-200H-7130 

5-2-88 

A2 

532  nm  Upwelling  Light, 
Biospherical  Instr. 

MCP-200H-7131 

5-2-88 

A3 

Transmissiometer, 

Sea  Tech,  Inc. 

165 

12-31-86 

A4 

Fluorometer,  Sea  Mar 

Tech,  Inc.  (signal) 

6000AR-235 

None 

A5 

Fluorometer,  Sea  Mar 

Tech,  Inc.  (scale) 

6000AR-235 

None 

A6 

488  nm  Downwelling 

Light,  Biospherical  Instr. 

QCP-200LM-7101 

2-22-86 

A7 

Echo  Sounder,  Ulvertech, 

LTD  . 

205 

unknown 

★ 

+ 


TUPS  on-board  computer  channel 
Date  of  last  calibration  if  known 


TESS  also  contained  two  narrow-band  irradiance  sensors  with 
fixed  wavelengths  of  441  and  488  nm.  Table  2  lists  the 


sensors  used  during  the  Panama  City  experiment . 

Table  2.  TESS  Sensor  Configuration,  June  1988 


Channel 

Sensor 

4 

Voltage  Reference,  2.5V,  Analog  Devices  AD580M 

5 

Pyroheliometer  1,  light  bulb,  Eppley 
Instruments,  Model  50,  S/N  3039 

8 

Pyroheliometer  2,  hemisphere,  Eppley 
Instruments,  Model  PSP,  S/N  8022D1 

11 

488  nm  Light  (linear),  Biospherical 

Instruments,  Model  QCP-200HM-488,  S/N  7105 

12 

441  nm  Light  (log),  Biospherical  Instrument, 
Model  QCP-200LM-441,  S/N  7102 

2.5  TUPS  At-Sea  Data  Collection 

All  sensors  from  the  TUPS  tow  vehicle  output  their 
signals  to  an  onboard  computer.  The  Sea  Bird  temperature  and 
conductivity  sensors  produce  a  frequency  output  and  all  other 
sensors  have  analog  (voltage)  outputs.  These  signals  are 
converted  by  the  TUPS  computer  to  a  hexadecimal  code  which  is 
sent  in  ASCII  format  to  the  computer  aboard  the  vessel.  The 
shipboard  computer  used  to  integrate  and  collect  data  from 
the  TUPS  computer  was  a  Digital  Equipment  Corporation 
Professional  350  (DEC  PRO-350)  running  under  the  P/OS 
operating  system.  The  DEC  PRO-350  collected  a  line  of  data 


from  the  TUPS  computer  every  5-6  seconds,  added  a  date  and 
time  to  the  data  line  and  wrote  the  data  line  to  the 
computer's  hard  disk.  This  operation  was  controlled  by  a 
BASIC  program  (SLOGTUFS)  which  collects  the  data,  calculates 
and  displays  the  results  in  engineering  units  and  writes  the 
raw  data  to  disk. 

2.6  TESS  At-Sea  Data  Collection 

The  data  collected  from  the  TUPS  environmental  sensor 
system  (TESS)  include  navigation  information  as  well  as 
light  data.  TESS  data  are  collected  with  a  Zenith  Data 
Systems  Model  121  (Z-121)  computer  running  MS/DOS  version 
2.17.  The  navigation  information  is  transmitted  to  the  Z-121 
computer  via  an  RS-232  interface.  The  navigation  data  is 
collected  from  an  INTERNAV  LC-300  LORAN-C  which  outputs 
position,  time  delays  and  calculated  speed  and  heading  every 
12  seconds.  The  LORAN  output  is  used  as  a  trigger  for  the  Z- 
121  BASIC  data  logging  program  (LOGTESS)  .  When  the 
navigation  data  is  sent  to  the  Z-121,  the  computer  records 
the  day  and  time  and  samples  the  A/D  board  that  is  receiving 
data  from  the  TESS  light  sensors.  Each  light  sensor  is  read 
five  times  and  an  average  is  determined.  The  time  and 
position  information  along  with  the  average  light  voltages 
are  displayed  on  the  Z-121  screen  and  simultaneously  written 
to  the  floppy  disk  in  ASCII  format. 
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2.7  Data  Collection  Problems 

1  Relatively  few  problems  occurred  during  collection  of 

the  data.  The  LORAN-C  positions  off  Panama  City,  Florida, 
are  very  good  due  to  the  near-perpendicular  crossings  of  the 

►  time  delay  lines.  Most  of  the  TUPS  sensors  were  well- 
calibrated  and  produced  acceptable  data  with  a  few  exceptions 
(e.g.  the  fluorometer)  that  could  be  connected  by  judicious 

*  filtering  of  the  data. 

The  one  major  problem  did  not  involve  the  TUPS  or  TESS 
equipment,  but  rather  the  loss  of  ship's  power  at  1615Z  on 

►  June  20,  1988.  This  event  is  described  in  section  2.2.  The 

loss  of  power  caused  the  computer  disk  files  on  the  DEC-350 
(hard  disk)  and  the  Z-121  (floppy  disk)  not  to  be  closed 

)  properly.  Although  the  data  had  been  written  to  the  disks 

immediately  after  it  was  collected,  it  could  not  be  read  in 
the  normal  fashion  because  location  blocks  for  the  file  had 
)  not  been  written  to  the  disk.  These  data  were  later 

recovered  using  a  technique  that  is  described  in  section  3.3. 

►  3.0  DATA  ANALYSIS 
3.1  Inspection  of  Data  Collected 

Immediately  after  the  Panama  City  experiment,  the  TUPS 
^  and  TESS  data  files  were  transferred  from  the  disk  on  which 

they  were  recorded  (DEC  350  or  Z-121)  to  the  GERG  VAX  for 
processing.  The  files  were  scanned  to  make  sure  each  data 
^  line  was  complete  and  that  each  file  had  ended  with  a 

complete  line  of  data.  Incomplete  data  lines  were  either 
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removed  or  edited  if  the  error  was  obvious.  Most  of  the  data 


files  were  in  excellent  condition.  It  was  obvious  however, 
that  there  was  a  problem  with  spiking  in  the  fluorometer  data 
that  would  need  attending  to. 

3.2  Problems  with  Collected  Data 

After  the  initial  inspection,  the  TUPS  and  TESS  files 
were  converted  to  engineering  units .  The  programs  used  are 
described  in  section  3.5.  A  quick  look  at  the  raw  data 
indicated  that  there  were  significant  problems  with  the  data 
from  downwelling  light  sensor  in  TUPS  and  the  488  nm  TESS 
(atmospheric)  sensor.  The  downwelling  light  sensor  (looking 
up)  in  TUPS  was  supposed  to  be  a  Biospherical  Instruments 
400-700  nm  broadband  sensor  that  measures  phot csynthet ical ly- 
active  radiation  (PAR) .  Inadvertently,  a  narrow  band  488  nm 
sensor  (Biospherical  Instruments,  Inc.  Model  QCP-200LM-7101) 
was  installed  in  its  place.  This  sensor  was  intended  for  the 
TESS  system.  The  488  nm  sensor  in  TUPS  was  a  logarithmic 
sensor  intended  to  operate  in  air.  The  calibration  for  this 
unit  is  an  air  calibration  thus  the  values  reported  for  this 
sensor  are  questionable. 

An  incorrect  sensor  was  also  installed  in  the  TESS. 
Instead  of  installing  the  logarithmic  448  nm  sensor  in  TESS, 
a  high  gain  488  nm  sensor  (Biospherical  Instruments,  Inc. 
Model  QCP-200HM-488-7105)  had  been  installed  previously 
during  the  Spring  of  1988.  The  calibration  for  this  unit  was 
for  use  in  water  only,  thus  the  data  from  this  sensor  is  also 


questionable.  In  fact,  examination  of  the  time  series  plots 
indicates  that  the  sensor  was  saturated  almost  all  the  time 
during  this  three  day  experiment . 

A  more  significant  problem  was  the  fluorometer  data. 
About  half  of  the  data  points  were  unusable.  The  SeaMarTech 
fluorometer  can  produce  a  noisy  signal  due  to  the  capacitor 
discharge  when  its  strobe  light  is  flashing.  If  the  TUPS 
computer  samples  at  this  time,  an  errant  data  point  will  be 
recorded.  We  have  recommended  that  NORDA  add  an  electronic 
filter  to  their  fluorometer  to  correct  this  problem  in  the 
future.  There  is  also  a  problem  with  the  fluorometer  when  it 
is  changing  scales.  The  fluorometer  may  jump  between  scales 
rapidly  before  it  settles  on  one  scale.  When  the  TUPS 
computer  samples  during  this  period  of  scale  change  an  errant 
data  point  is  recorded.  Since  the  Sea  Mar  Tech  fluorometer 
does  not  change  scale  very  smoothly  this  process  also  added 
to  the  bad  fluorometer  data  recorded. 

There  were  only  a  few  problems  with  the  other  TUPS 
sensors.  When  the  vessel  ran  aground,  much  sediment  was 
stirred-up  into  the  water  column.  Sensors  that  work  by 
having  water  pass  through  them  (conductivity  and 
transmission)  are  affected  by  this  extra  suspended  sediment. 
Both  these  sensors  had  spikes  in  the  data  lines  when  the  ship 
ran  aground.  These  spikes  were  filtered  from  the  data. 


3.3  Recovery  of  Lost  data 

On  20  June  1988  the  generator  overheated  on  the  Captain 
Graydon  York  causing  the  unit  to  shut  down.  The  resultant 
power  failure  caused  TUPS  data  files,  written  on  the  DEC  PRO- 
350  hard  disk,  not  to  be  closed.  Although  the  data  were 
written  on  the  disk,  the  directory  entry  for  the  data  files 
was  not  complete  and  the  files  could  not  be  read  in  the  usual 
fashion.  Features  for  the  recovery  of  the  data  were  not 
available  in  the  P/OS  operating  system  of  the  DEC  PRO-350 
computer.  We  undertook  the  job  of  recovering  these  lost  files 
so  the  TUPS  data  from  June  20,  1988  could  be  saved. 

Our  strategy  for  recovering  those  files  was  as  follows: 

a.  Boot  an  RT-11  system  from  a  floppy  disk  on  the  PRO- 
350.  RT-11  is  another  DEC  operating  system  that  can  be  used 
on  the  PRO-350.  FORTRAN  subroutines  are  available  in  the  RT- 
11  system  library  that  can  bypass  normal  directory  usage  and 
allow  individual  blocks  of  data  to  be  read  anywhere  on  the 
hard  disk. 

b.  On  another  DEC  LSI-11  computer  available  at;  TAMU- 
GERG,  develop  a  FORTRAN  program  that  could  be  operated  under 
RT-11  and  could  search  the  disk  and  recover  blocks  containing 
Panama  City  data.  All  of  the  lost  TUPS  data  contained  the 
string  "20- JUN-88"  so  the  program  was  instructed  to  copy  all 
blocks  that  contained  this  string  to  a  new  file  on  a  RT-11 
floppy  disk.  The  program,  named  READ350,  is  listed  in 
Appendix  A. 
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c.  The  program  READ350  was  run  on  the  PRO-350  and  all 
data  recovered  to  a  RT-11  floppy  disk. 

d.  The  data  on  the  RT-11  floppy  disk  was  transferred  to 
the  VAX  using  the  DEC  LSI-11  computer  available  at  TAMU-GERG. 

e.  The  data  was  edited  to  clean  up  a  few  bad  characters. 
Essentially  all  lost  data  was  recovered  in  this  manner. 

The  TESS  data  collected  on  the  Zenith  121  computer  was 
recovered  by  Robert  A.  Arnone  at  NORDA  using  a  "Brown  Bag" 
data  recovery  program.  This  recovery  process  provided 
multiple  ASCII  files  of  the  data  which  were  slightly 
overlapping.  The  recovered  data  files  were  transferred  to 
the  GERG  VAX  merged  into  a  single  file  and  edited  to  remove 
all  duplicate  data  records  produced  by  the  recovery  program. 
This  recovery  of  the  TESS  data  from  June  20,  1988  provided  a 
clean  TESS  file  for  processing  along  with  the  files  from  the 
other  two  days. 

3.4  Data  Analysis  Protocol 

TESS  data  records  contain  time  and  voltages  in  ASCII. 
TUPS  data  records  contain  time  and  date  in  ASCII  together 
with  voltages  and  frequencies  in  a  hexadecimal  format.  Some 
sensors  are  inherently  noisy.  TESS  and  TUPS  data  are  not 
collected  at  exactly  the  same  time.  The  goal  of  the  data 
analysis  protocol  is  to  convert  the  raw  TESS  and  TUPS  data 
into  a  unified  data  set  containing  time  series  of  all 
measured  variables  in  engineering  units  with  erroneous  data 
and  noise  removed  as  much  as  possible. 
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Analysis  proceeds  in  a  number  of  steps  with  intermediate 
results  stored  as  files.  The  process  is  shown  schematically 
in  Figure  3  and  Figure  4.  Square  boxes  indicate  data  sets 
and  parallelograms  represent  computer  programs. 

First,  the  TESS  data  are  read  and  converted  into 
engineering  units  by  program  READJUN.  Formulas  used  to 
convert  voltages  to  engineering  units  are  best  found  in  the 
programs  in  Appendix  A. 

Second,  TUPS  data  are  converted  to  engineering  units  by 
program  CONVERT 8 8 . 

Third,  certain  noisy  measurements  (depth,  fluorescence) 
are  filtered  using  program  PRETABLE.  Rather  than  just 
applying  a  low  pass  filter,  PRETABLE  tries  to  make 
intelligent  decisions  about  what  data  points  are  noise  and 
removes  them  from  the  time  series. 

Fourth,  TESS  data  are  interpolated  to  the  time  of  the 
TUPS  data  and  interpolated  TESS  data  and  the  TUPS  data  are 
merged  using  program  MERGETT.  This  is  done  for  each  of  the 
three  days.  The  result  is  3  files  containing  the  desired 
smooth  data  sets  of  all  data  collected  on  each  day  of  the 
Panama  City  experiment. 

Fifth,  a  quick  look  data  plotting  program  (TTPLOT)  was 
written  to  read  the  merged  data  set  and  plot  each  time  series 
on  on  large  poster-sized  plot. 

Sixth,  a  SAS  program  (READ.SAS)  was  written  to  read  all 
final  data  sets  and  consolidate  them  into  one  SAS  data  set. 
SAS  data  sets  are  self  documenting  and  facilitate  the  use  of 


Figure  3.  Flow  diagram  of  data  processing  protocol  showing 
FORTRAN  programs  used  and  data  file  names. 


Figure  4.  Flow  diagram  of  data  processing  protocol  showing  SAS 
programs  used  and  data  file  names. 


powerful  SAS  statistical  routines.  Another  SAS  program  was 
written  to  extract  subsets  of  environmentally  significant 
variables  for  further  analysis. 

Seventh,  a  contouring  program  (CONTOUR)  was  used  to 
prepare  contour  plots  for  each  of  the  three  days  for  all 
environmentally  significant  variables. 

The  smooth  data  files  that  we  recommend  be  used  for 
further  processing  are  the  19JUN88.FIN,  20JUN88.FIN  and 
21 JUN88 . FIN  files  produced  by  the  MERGETT  program.  Table  3 
gives  a  listing  of  the  position  in  the  record  (array  number) 
of  each  of  the  data  items,  the  descriptive  name  of  the  data 
item,  the  SAS  name  used  for  the  correlation  analysis  and  the 
units  of  each  data  item  as  reported. 

3.5  Description  of  Computer  Programs 

The  following  is  a  brief  description  of  the  programs 
that  were  used  to  analyze  the  data  and  to  prepare  the  figures 
in  this  report.  The  listings  of  the  programs  are  presented 
in  Appendix  A  for  information  purposes  only.  The  .listing 
does  not  imply  any  transfer  of  ownership  or  guarantee  that 
the  program  will  operate  on  any  data  sets  or  computers  other 
than  those  actually  used  in  the  course  of  this  work.  Some  of 
these  programs  are  based  on  previously-developed  programs. 
They  use  subroutines  of  commercial  or  public  domain  origin 
available  on  the  TAMU-GERG  VAX. 
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Table  3.  Position  in  data  record  (from  XXJUN88.FIN  files), 

descriptive  name,  SAS  name  and  units  for  data  presented 
in  this  report  and  on  the  data  tape  supplied  to  NORDA. 


Position  Descriptive  Name 


SAS  Name  Units 


» 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 


Year  -  - 

Month  -  - 

Day  -  - 

Time  HH:MM:SS(UT) 

Seconds  Past  Midnight  sec 

Latitude  deg 

Longitude  deg 

Time  Delay  1  flsec 

Time  Delay  2  Usee 

Heading  deg 

Speed  -  '  knots 

Voltage  Reference  volts 

441  nm  light  (TESS)  TESS441  )lW/cm2sec  nm 

488  nm  light  (TESS)  TESS488  HW/cm2  nm 

Pyroheliometer  1  TESSPYR1  |XW/cm2sec*l0-4 

Pyroheliometer  2  TESSPYR2  Ji.W/cm2sec*10~4 

Temperature  TEMPER  °C 

Salinity  SALIN  PSU 

%  Transmission  PTRANS  % 

Fluorescence  FLUOR  rel.  units 

4  65  nm  Upwelling  Light  L4  65NM  jlW/CM2/nm 

507  nm  Upwelling  Light  L507NM  ^W/CM2/nm 

532  nm  Upwelling  Light  L532NM  fj.W/CM2/nm 

Transmissometer  Voltage  Volts 

Fluorometer  Signal  Voltage  -  Volts 

Fluorometer  Scale  Voltage  -  Volts 

488  nm  Downselling  Light  L488NM  p.W/cm2sec  nm 

Depth  DEPTH  Meters 
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1.  CONVERT88  (convert  raw  TUPS  files) . 

This  program  reads  the  raw  hexadecimal  data  file 
produced  by  the  TUPS  as  configured  during  June  1988  and 
converts  data  to  engineering  units.  The  output  files 
produced  by  this  program  have  the  extension  "CNV" . 

2.  READJUN  (read  June  {TESS}  files). 

This  program  reads  the  files  produced  by  the  TESS 
program  (as  configured  during  June  1988)  and  produces  a  clean 
data  file  in  engineering  units.  The  output  files  produced  by 
this  program  have  the  extension  "RDJ" . 

3.  PRETABLE  (preprocess  TUPS  file). 

This  program  reads  the  TUPS  files  produced  by  CONVERT88 
and  filters  the  data.  A  special  filter  (which  does  not 
average-in  fliers)  is  used  for  the  fluorescence,  depth,  and 
some  other  data.  The  output  files  produced  by  this  program 
have  the  extension  "PTO” . 

4.  MERGETT  (merge  TESS  and  TUPS  data). 

This  program  reads  the  output  from  READJUN  (TESS)  and 
PRETABLE  (TUPS)  and  produces  a  uniform  data  set  containing 
TESS  and  TUPS  data.  TESS  data  is  interpolated  to  TUPS  times. 
The  output  files  produced  by  this  program  have  the  extension 
"FIN" . 

5.  TTPLOT  (TESS  TUPS  plot). 

This  program  uses  NCAR  subroutines  to  plot  all  elements 
of  the  unified  data  set  against  time.  A  large  format  plot 
showing  all  data  is  output  on  32  in  by  42  in  inch  paper. 
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6.  CONTOUR  (contour  13  TESS-TUPS  variables) . 

This  program  contours  data  along  the  cruise  track.  The 
cruise  track  is  show  as  a  dotted  line.  Land  areas  are  also 
displayed.  A  masking  algorithm  is  used  to  produce  contours 
only  within  a  certain  distance  of  the  cruise  track.  The  same 
algorithm  is  used  to  mask  out  contours  that  might  be  drawn 
over  land. 

7 .  Ml  (map  one)  . 

This  program  is  used  to  convert  digitizer  (x,y)  data  to 
latitude-longitude  coordinates.  A  latitude-longitude  grid  is 
first  digitized  from  20  known  positions  to  produce  a 
calibration  field.  Then  the  data  locations  are  digitized.  The 
program  does  a  third-order  orthogonal  polynomial  fit  to  the 
calibration  points  and  then  uses  the  resultant  coefficients 
to  calculate  latitude  and  longitude  for  the  data  locations. 
The  program  was  used  to  provide  the  digitized  island  shown  in 
the  contour  maps . 

8.  READ . SAS  (read  all  final  data  sets) 

This  program  reads  final,  merged  data  sets  (date. fin) 
and  converts  all  Panama  City  data  into  a  SAS  data  set. 

9.  PUTALL .  SAS  (extract  subset  of  environmental  data  for 
plotting  to  ASCII  file) 

This  is  a  SAS  program  that  can  use  simple  statements  to 
make  an  ASCII  file  to  be  used  by  the  plotting  programs. 

10.  CORR.SAS  (make  correlation  matrix) 

This  program  produces  a  correlation  matrix  of  all  the 


Panama  City  data. 


11.  READ350  (read  and  recover  PRO-350  data) 

This  is  an  RT-11  program  that  uses  RT-11  system 
subroutine  calls  to  read  the  hard  disk  on  the  PRO-350  and 
recover  blocks  of  lost  data.  These  data  could  not  be  read  in 
the  usual  fashion  because  their  directly  entry  was  incomplete 
due  to  a  power  failure  while  the  data  were  being  written  to 
the  PRO-350  disk. 

12.  PLOT . S AS  (make  x  vs .  y  plot) 

This  SAS  program  plots  the  inverse  of  NUR  against  depth. 
Inverse  NUR  (INORM_UP  =  TESSP YR1/L507NM)  is  plotted  against 
depth  for  each  day  for  all  depths  less  than  12  meters 

4.0  DATA  PRESENTATION 
4.1  Time  Series  Plots 

Time  series  data  for  each  day  were  plotted  on  42"  by  32" 
paper  to  facilitate  examination  of  the  data.  Page  size 
versions  of  these  plots  are  shown  in  Figure  5,  6  and  7  for 
June  19,  June  20,  and  June  21.  Each  time  series  is  scaled 
such  that  it  fills  its  own  box.  The  minimum  and  maximum 
value  fox  each  box  is  indicated  near  the  right  hand  edge. 
Time  starts  on  even  hours.  For  example.  Figure  6  shows  1500Z 
tc  1700Z.  It  is  apparent  that  the  TUPS  upwelling  irradiance 
sensors  (Panels  7,  8,  and  9  from  the  top)  clearly  measure  the 
increase  in  upwelling  irradiance  in  shallow  water.  Also 
evident  is  the  general  inverse  relationship  between 
fluorescence  and  transmission. 
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4.2  Contour  Maps  for  Each  Day 

Appendix  B  contains  contour  plots  for  13  environmental 
variables  for  each  day  of  the  cruise.  These  are  intended  to 
give  a  quick  look  at  the  spatial  quality  of  the  data.  Where 
cruise  traces  are  relatively  uniformly  spaced  and  there  are 
more  than  one  track  as  on  June  19,  the  data  such  as  depth 
contour  nicely.  Passage  of  clouds  show  up  dramatically  as 
holes  in  the  light  field. 

4.3  SAS  Analysis  of  Total  Data  Set 

As  mentioned  in  sections  3.4  and  3.5,  all  data  were 
placed  in  a  SAS  data  file  to  facilitate  statistical  analysis. 
The  Pearson  correlation  matrix  was  calculated  for  the  entire 
data  set  (all  3  days)  for  the  flowing  variables: 

•  TESS  light  sensors:  441  nm,  448  nm,  Pyrometer  1, 
Pyrometer  2 

•  TUPS  sensors:  Temperature,  Salinity,  Depth, 
Fluorescence,  Percent  Transmission 

•  TUPS  light  sensors:  465  nm,  488  nm,  507  nm,  532  nm. 

This  matrix  is  shown  in  Table  4.  Correlation  coefficients 
are  highlighted  when  their  absolute  value  exceeds  0.60  to 
draw  attention  to  the  most  significant  correlations. 

Three  of  the  four  TESS  light  sensors  are  highly 
correlated  with  correlation  coefficients  exceeding  0.97.  The 
488  nm  TESS  sensor  has  a  low  correlation  with  the  other  three 
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TESS  light  sensors.  This  provides  additional  evidence  that 
the  TESS  488  nm  sensor  was  not  working  properly. 

Fluorescence  is  negatively  correlated  with  transmission 
which  implies  that  the  principle  light  scatterers  in  the  area 
are  phytoplankton.  Fluorescence  is  negatively  correlated 
with  the  three  light  sensors  which  probably  reflects  the 
increasing  efficiency  of  plankton  fluorescence  with 
decreasing  ambient  light.  Fluorescence  is  negatively 
correlated  with  salinity  which  implies  the  fresher  bay  waters 
contain  more  chlorophyll. 

Salinity  and  temperature  are  negatively  correlated;  the 
bay  waters  are  warmer  and  fresher  than  the  open  ocean  waters 
off  the  beach.  Three  TUPS  light  sensors  measuring  upwelling 
light  are  highly  correlated. 

We  made  a  quick  look  assessment  of  the  depth  and  light 
data  to  see  how  well  water  depth  could  be  predicted  by  the 
intensity  of  upwelling  irradiation.  First  we  normalized  the 
upwelling  intensity  to  the  downward  radiation  measured  by  the 
on  board  TESS  sensors.  Since  the  three  TESS  light 
measurements  were  highly  intercorrelated  and  since  the  three 
TUPS  measurements  of  upwelling  irradiation  were  also  highly 
correlated,  we  felt  justified  in  choosing  one  of  each  for 
this  normalization.  We  defined  normalized  upwelling 
radiation,  NUR,  to  be  the  TUPS  507  nm  sensor  data  divided  by 
the  TESS  Pyroheliometer  1. 

The  inverse  of  NUR  is  plotted  against  depth  in  meters  in 
Figure  8.  This  is  a  SAS  plot  where  the  letters  indicate  the 
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number  of  data  items  under  each  letter.  An  A  represents  one 
data  item,  a  B  represents  two  data  items,  and  so  forth.  A  Z 
indicates  26  or  more  data  items.  Figure  8  shows  that  for  all 
the  data  collected  on  June  19,  depth  is  roughly  proportional 
to  the  inverse  of  NUR.  Figure  9  show  the  same  for  the  June  20 
data.  Again  depth  seems  predictable  by  the  inverse  of  the 
NUR. 

Figure  10  for  June  21  show  a  very  different  picture  for 
the  data  collected  behind  the  barrier  islands.  No  simple 
relationship  exists  between  NUR  and  water  depth.  Bottom  types 
encountered  on  June  21  included  mud,  bright  sand  and  grass. 
Bottom  types  encountered  on  June  19  and  20  were  mostly  bright 
sand.  Bottom  reflectance  may  play  a  role  in  distorting  the 
relationship  between  NUR  and  water  depth.  Other  environmental 
variables  such  as  turbidity  and  pigment  concentrations  most 
certainly  also  play  a  role.  The  data,  of  course,  warrant 
further  analysis. 

4.4  Description  of  Deliverables 

Provided  to  NORDA  in  addition  to  this  report  were  three 
42"  x  32"  versions  of  Figures  5,  6  and  7.  Also  provided  was  a 
nine  track  magnetic  tape  in  VAX  BACKUP  format  that  contains 
all  data  files  described  in  this  report,  plus  copies  of  all 
FORTRAN  and  SAS  programs.  Data  files  generally  are  ASCII 
FORTRAN  output  files  with  a  carriage  return  line  feed  pair 
(CR-LF)  at  the  end  of  each  line.  The  BACKUP  log  that  contains 
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Figure  8.  Plot  of  inverse  of  normalized  upwelling  irradiance  versus  depth  on  19  June  1988. 
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Figure  10.  Plot  of  inverse  of  normalized  upwelling  irradiance  versus  depth  on  21  June  1988. 


a  listing  of  the  names  of  all  files  on  the  tape  is  provided 
as  Appendix  C. 


5.0  CONCLUSIONS 

All  TUPS  and  TESS  data  collected  from  the  NORDA  Panama 
City  experiment  has  been  processed,  merged,  and  placed  in  a 
form  amenable  for  further  analysis.  The  data  set  is 
relatively  complete,  as  we-  were  able  to  recover  the  data 
"lost"  due  to  a  ship  power  failure  on  June  20,  1988.  From  the 
combined  TUPS/TESS  data,  we  produced  large-scale  time  series 
plots  showing  all  variables.  These  plots  are  excellent  tools 
for  examining  the  quality  of  the  data  set  and  for  taking  a 
quick  look  at  parameter-parameter  relationships.  All  light 
data  appear  to  be  of  high  quality,  with  the  possible 
exception  of  the  488  nm  downwelling  light  data.  The  TUPS 
depth  sensor  seemed  to  work  exceptionally  well. 

The  TUPS/TESS  system  appears  to  be  quite  useful  for 
examining  the  optical-depth  relationship  in  shallow  water. 
Water  depth  can  be  inferred  from  measurements  of  upwelling 
irradiance  under  some  conditions.  Preliminary  examination  of 
this  data  set  indicates  the  relationship  holds  well  in  areas 
of  bright  sand  (e.g.  south  of  Shell  Island),  but  the  relation 
of  light  to  depth  is  more  complicated  in  variable  bottom  and 
water  clarity  conditions.  The  broad  suite  of  light  and  water 
quality  sensors  in  TUPS  makes  it  ideal  for  resolving  such 
problems . 
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COMPUTER  PROGRAM  LISTINGS 


ft  ft  ft  ft  ft  ft  ft  ft  ft~ft  c 
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OflRAM  C0NVRT88 

PURjPjOSE  1)  reads  in  raw  TUPS  data  files  in  hex 

2)  converts  data  to  engr.  units  with  proper  calib. 
TjHjOR  D.A.UIESENBURG 

MjE|  2/9/88 


REVi 


€ 


*t^****  **********  ****************************************************l* 

CHARACTER* 3  MONTH 
CHARACTER*8  CHTIME.HHMMSS 
CHARACTER*9  CHDATE , DDMMYY 
CHARACTER*30  FILEN 
CHARACTER*34  CNVDAT , RAVDAT 
CHARACTER*60  RRAW 
CHARACTER*? 6  BUFFER 
CHARACTER*4  ANALG(8) 

CHARACTER*6  FRQ(8) 


SION 

INPUT  FILE  IS  GIVEN  AND  OUTPUT  FILE  OF  SAME  NAME . CNV  PRODUCED 


NORMAN  GUINASSO 
DEC.  88 


REAL* 4  V(8) 

REAL*8  TIME1( 10000), TIMEA( 10000) 


P 


INTEGER*4  STR$FIND_FIRST_IN_SET 
COMMON  /DAYS/  IDAY,IMON,IYR, jhr ,j rain 
COMMON  /DATIM/  CHDATE, CHTIME 
COMMON  /  RAW  /  ANLG , FRQ 
COMMON  /  DAY  /  DMY,HMS 

COMMON  /  CNVRT  /  COND,SAL78, TEMP, VOLT, FREQ, XL4,XL8,PSI, 
DEPTH , TRANS , TLIGHT , CHLA , MINUTE , RTIME 


oldv5=1.0  !  initialize  value  used  for  fluor.  sig.  cleaning 


P 


c 


ioq 


W 


VRITE(6, 100) 

FORMAT ('  ENTER  FILENAME  FOR  CURRENT  RUN',/, 

'  INCLUDING  EXTENSION  (<=  30  CHARS)',/) 
READ(*, ' (a) ' )  FILEN 
rawdat=f ilen 


call  chg_ext(rawdat ,cnvdat , 'UNF' ) 

OPEN(UNIT=18, FILE=RAWDAT, STATUS=' OLD' ) 

OPEN(UNIT=19 , FI LE=CNVDAT , STATUS= ' NEW ' , form=' unformat  ted ' , DISP= ' DELETE' ) 
call  chg  ext(rawdat,cnvdat, 'CNV' ) 


i=0 

do  while( . true. ) 

read(18, ' (a)' , iostat=ios)  buffer 
if(ios.eq.-l)  goto  200 
i  =  i  +  l 

J=STR$f ind_f irst_in_set(buf fer, ' : ' ) 
read(buffer( j-2: j+2) , ' (i2, lx, i2) ' )  ih, im 
timel(i)  =  60.d0*(60.d0*ih  +  im) 
jhr  =ih 

jmin  =im 

the  seconds  are  all  00  on  TUPE  data  files 


Page  2 


I 


I 


> 
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j=str$f ind_first_in_se t( buffer , '  -'  ) 

read(buffer( j-2: j+6) , ' ( i 2 , lx,a3, lx, i2) ' )  iday, month, iyr 

IMON  =  IDECMONTH( MONTH) 

inum=i 

RRAW  =  BUFFER( 1:60) 

CALL  CONVRT ( RRAW , i ) 

enddo 


1200 


call  timefix(timel, timea, inum) 
rewind( 19) 

CALL  DATE(DDMMYY) 

CALL  TIME(HHMMSS) 
open( 20, name=cnvdat , s tatus= ' new' , 
RECL=255, carriagecontrol=' list' ) 


► 


> 


» 


I 


I 


I 


C 

c 


i5d 


i6d 


1 

2 


3 

4 


WRITE(20, 150)  CNVDAT , DDMM YY , HHMMSS 

FORMAT ( '  FILE:  ',A35,5X,'  CREATED:  ' ,A9,2X,A8,/) 

WRITE(20, 160) 

FORMAT ( IX , 4HDATE , 5X , 5HSTIME , 4X , 5HZTIME , 5X , 4HTEMP , 3X , 5HSAL78, 

4X , 5HCTRAN , 3X , 6HCHL-FL , 5X , 5H465NM , 3X , 5H507NM , 3X , 5H532NM , 

3X, 6HTRANSV, 2X, 6HFLSIGV, 2X, 6HFLSCLV, 2X, 6H488-UP , 3X , 5HDEPTH , 
2X,7('-'),3X,8('-'),1X,5('-'),2X,5('-'),2X,6('-'), 
2X,5('-'),2X,6('-'),2X,5('-'),2X,5('-')) 


/> 


* 

★ 

* 


DO  WHILE  (.TRUE.) 

read(19,end=500) 

i, IDAY, IMON, IYR, jhr.jmin, TEMP, SAL78, TRANS, CHLA, (V(K) ,K=1 
write(20, 1000,IOSTAT=IOS)  iyr , imon, iday , timea(i) , 
jhr.jmin, 

TEMP, SAL78, TRANS, CHLA, (V(K),K=1, 8) 


,0) 


ENDDO 


500 

1000 


H 


close(20) 

CLOSE(UNIT=18) 

CLOSE(UNIT=19) 

FORMAT(I2.2, ' / ' , 12.2, f / ' ,12. 2, lx, f8. 1 , 2x, i2 . 2, ' : ' , i2. 2 , ' :00' , lx, 
F6.2,F8.3,F8.2,F9.2,1X, ' | ' ,6F8.3,F9.3,1X,F7.2) 


STOP 

END 


W 


integer*4  function  idecmonth(chr) 
character  *(*)  chr 
character*3  months(12) 
data  months  /  ' JAN' , 'FEB' , 'MAR' , 'APR' 

'JUL',' AUG', 'SEP', 'OCT' 
call  str$upcase(chr,chr) 

DO  1=1,12 

IF(CHR.EQ.MONTHS(I))  THEN 
IDECMONTH=I 
RETURN 


'MAY' 

'NOV' 


'JUN'  , 
'DEC'/ 


ENDIF 


ENDDO 

TYPE  *,'  IDECMONTH-W-BAD  MONTH' 

RETURN 

END 


ic 


SUBROUTINE  CONVRT ( RRAW, i ) 

CNVRT  opens  existing  TUPS'  data  files,  one  at  a  time,  and  first 
divides  the  strings  into  frequencies  and  analogs,  then  converts 
the  strings  to  various  data  parameters,  which  are  written  to 
another  output  file. 

modified  at  Texas  A&M  Univeristy  by  Guinasso  and  Wiesenburg 


I 


cvnrn  rtTi  rt  <n  <n  n  rrrrrr T5r 
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June  1988  —  for  Panama  City  Cruise 

CALLED  FROM  COLLCT  AND  CALLS  CNVTLT  AND  PARSE 

CHARACTER*60  RRAV 
CHARACTER*4  ANLG(8) 

CHARACTER*6  FRQ(8) 

CHARACTER*8  CHTIME 
CHARACTER*9  CHDATE 
REAL  V(8) 

i 

COMMON  /DAYS/  IDAY, IMON, IYR, jhr, jmin 
COMMON  /  DATIM/  CHDATE, CHTIME 

COMMON  /  CNVRT  /  COND,SAL78, TEMP, VOLT, FREQ, XL4,XL8,PSI, 
1  DEPTH , TRANS , TLIGHT , CHLA , MINUTE , RTIME 

COMMON  /  RAW  /  ANLG, FRQ 

DATA  XLFACT  /0.121E17/ 

Blank  out  the  ANLG  and  FRQ  arrays 

CALL  ZERO 

Divides  TUPS  string  into  frequencies  and  analogs 
CALL  DIVIDE(RRAV) 


Convert  data  parameters 


P 


CALL  CNVTMP 
CALL  CNVCND 
CALL  CNVPSI 


!  temperature  -  1st  frequency 
! conductivity  -  2nd  frequency 
!  pressure  -  3rd  frequency 


P 


FORWARD  LIGHT  SENSOR  IN  FIRST  CHANNEL  0-5  VOLTS 
CALL  CNWLT(ANLG(1)  ) 

V(l)  m  VOLT*2.0833  !  Factor  for  2.4  volt/5  volt  conv. 

V(l)  .  V(l)*10.9459  !  CAL.  FACTOR  FOR  465  LIGHT  SENSOR  6/88 

!  UNITS  ARE  uV/cm2/nm 


P 


MIDDLE  LIGHT  SENSOR  IN 
CALL  CNWLT(ANLG(2)) 
V(2)  =  V0LT*2.0833 
V(2)  =  V(2)*10.3644 


SECOND  CHANNEL  0-5  VOLTS 

!  Factor  for  2.4  volt/5  volt  conv. 
!  CAL.  FACTOR  FOR  507  SENSOR  6/88 
!  UNITS  ARE  uV/cm2/nm 


P 


AFT  LIGHT  SENSOR  IN  THIRD  CHANNEL  0-5  VOLTS 
CALL  CNWLT(ANLG(3)) 

V(3)  =  V0LT*2.0833  !  Factor  for  2.4  volt/5  volt  conv. 

V(3)  =  V(3)*9.8568  !  CAL  FACTOR  FOR  532  LIGHT  SENSOR  6/88 

!  UNITS  ARE  uV/cm2/nm 


P 


TRANSMISSOMETER  IN  FOURTH  CHANNEL  0-5  VOLTS 
CALL  CNWLT(ANLG(4)  ) 

V(4)  =  V0LT*2.0833  !  Factor  for  2.4  volt/5  volt  conv. 

TRANS=100. *((4.738/4. 46 )*(V(4)-. 001 ))/5.  !  SENSOR  #165 


P 


FLUOROMETER  SIGNAL  IS  FIFTH  CHANNEL  0-10  VOLTS 
CALL  CNWLT(ANLG(5)) 

V(5)  =  V0LT*4.1667  !  Factor  for  2.4  volt/10  volt  conv. 


P 


check  for  bad  fluorometer  voltage  and  replace  if  bad 
if  (v(5)  .gt.  9.0)  then 
v(5)=oldv5 


U  *,  M  1  is 


i  /—DEC— 8  6  13  :26  :59  $OISKl : ( DENIS . PCITY . PGMS ] CONVRT88 . FOR; 2 

,  — —  _ *  _ _ _ _  _ _ _ • _ _ _ •  _ _ _ _ _ • _ _ _ _ _  •  _ _ _  • 

—  ,  —  •  —  •  —  •  — •  •  • - •  — •  — —  . - - -  •  a  —  —  —  —  J  —  - 

end  if 
oldv5=v(5) 

FLUOROMETER  SCALE  IS  SIXTH  CHANNEL  0-10  VOLTS 
CALL  CNWLT(ANLG(6)  ) 

V(6)  =  V0LT*4.1667  !  Factor  for  2.4  volt/10  volt  conv. 

CALCULATE  REL.  FLUORENCE  USING  SIGNAL  AND  SCALE  VALUES 
CHLA=V(5)*(2.0**ININT(V(6)))  !  Make  scale  volts  an  integer 

UP  LOOKING  LIGHT  SENSOR  IS  EIGHTH  CHANNEL  0-10  VOLTS 
UNITS  OF  uW/cm3/nm 
CALL  CNWLT(ANLG(7)) 

V(7)  =  VOLT*4.1667  !  Factor  for  2.4  volt/10  volt  conv. 

V(7)  =  EXP (V( 7) -0.900) *5. 88  !  448  LOG  SENSOR  USUALLY  IN  TESS 

****  THIS  IS  DIFFERENT  THAN  ALL  OTHER  TUPS  CRUISES  ***** 

ECHO  SOUNDER  IS  SEVENTH  CHANNEL  0-10  VOLTS 
CALL  CNVVLT(ANLG(8) ) 

V(8)  =  V0LT*4.1667  !  Factor  for  2.4  volt/10  volt  conv. 

V(8)  =  (V(8)*10.0)+1.0  !  meters=volts*10  +  1.0  m  tow  depth 

VRITE(19) 

i ,IDAY,IMON,IYR, j hr ,j min, TEMP, SAL78, TRANS, CHLA, (V(K) ,K=1,8) 

RETURN 

END 


SUBROUTINE  ZERO 

CHARACTER* 6  F(8) 
CHARACTER*4  A(8) 

COMMON  /  RAW  /  A,F 

DO  10  J=l,8 
F(J)  =  '  ' 

A(J)  =  '  ' 

10  CONTINUE 

RETURN 

END 


SUBROUTINE  HEXDEC(F) 

CHARACTER*4  0 
CHARACTER*6  F 
INTEGER*2  G 
CHARACTER *1  CM,NUM 

COMMON  /  CNVRT  /  COND, SAL78, TEMP, VOLT, FREQ, XL4,XL8,PSI, 
1  DEPTH, TRANS, TLIGHT, CHLA 

CM  =  F(lsl) 

XM  =  XMULT(CM) 

IF  (XM.NE.O. )  THEN 
DO  50  1=3,6 

0(I-2:I-2)  =  F(I:I) 

50  CONTINUE 
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iod 


c 

c*** 

c 


LENGTH  =  4 
XN  =  0. 

P  =  1. 

DO  100  I = LENGTH, 1,- 
G=ICHAR(0( 1:1)) 
IF(G.GT.64)  THEN 
D=FL0AT(G)-55. 
ELSE 

NUM  =  0(1:1) 

IF  (NUM.EQ. 'O' ) 


1 


D  =  0.  !  ADDED  BY  K.D. SAUNDERS  9/25/86 


D=FLOAT(G) 

END  IF 

XN=XN+D*P 

P=P*16. 

CONTINUE 

FREQ=10.0E06*XM/XN 

ELSE 

FREQ=0. 

END  IF 

[END  OF  SUBROUTINE  HEXDEC 

RETURN 
END 


IF 

(NUM.EQ. '1' 

) 

D 

= 

1. 

IF 

(NUM.EQ.  '2' 

) 

D 

= 

2. 

IF 

(NUM.EQ.  '3' 

) 

D 

= 

3. 

IF 

(NUM.EQ. '4' 

) 

D 

= 

4. 

IF 

(NUM.EQ.  '5' 

) 

D 

= 

5. 

IF 

(NUM.EQ. '6' 

) 

D 

= 

6. 

IF 

(NUM.EQ.  '7' 

) 

D 

= 

7. 

IF 

(NUM.EQ.  '8' 

) 

D 

= 

8. 

IF 

(NUM.EQ. '9' 

) 

D 

= 

9. 

FUNCTION  XMULT(CM) 

CHARACTER* 1  CM 

IF(CM.EQ.'O')  THEN 
XMULT=253. 

ELSE  IF(CM.EQ. ' 1' )  THEN 
XMULT=126. 

ELSE  IF(CM.EQ. ' 2' )  THEN 
XMULT=84. 

ELSE  IF(CM.EQ. '  3' )  THEN 
XMULT=63. 

ELSE  IF(CM.EQ. '4' )  THEN 
XMULT=50. 

ELSE  IF(CM.EQ. '5' )  THEN 
XMULT=42. 

ELSE  IF(CM.EQ. ' 6' )  THEN 
XMULT=36. 

ELSE  IF(CM.EQ. )  THEN 
XMULT=31. 

ELSE  IF(CM.EQ. '8' )  THEN 
XMULT=28. 

ELSE  IF(CM. EQ. ' 9' )  THEN 
XMULT=25. 

ELSE  IF(CM.EQ. 'A' )  THEN 
XMULT=23. 


O  U  O  O  O 
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C 

c*** 

c 


(ELSE  IF(CM.EQ.'B')  THEN 
XMULT=21 . 

(ELSE  IF(CM.EQ.'C')  THEN 
XMULT=19 . 

|ELSE  IF(CM.EQ.'D')  THEN 
XMULT=18. 

(ELSE  IF(CM.EQ.'E')  THEN 
XMULT=16 . 

pELSE  IF(CM.EQ. 'F' )  THEN 
XMULT=15. 

(ELSE 

XMULT=0. 

WRITE(6,*)  CM,'  *N0  MULTIPLIER  CALCULATED,  XMULT  SET  TO  ZERO' 
END  IF 

|END  OF  FUNCTION  XMULT 

RETURN 

END 


SUBROUTINE  CNVTMP 

CHARACTER*4  ANLG(8) 

CHARACTER*6  FRQ(8) 

COMMON  /  RAW  /  ANLG , FRQ 

COMMON  /  CNVRT  /  COND , SAL78 , TEMP , VOLT , FREQ , XL4 , XL8 , PSI , 

DEPTH , TRANS , TLIGHT , CHLA 

BIRD  CALIBRATION  DATA  FOR  SENSOR  #632 

A  /3.67517928E-03/,  B  /6.01320919E-04/,  „  !  SEABIRD  SENSOR 
C  / 1 . 59454020E-05/ ,  D  /2.59063611E-06/  !  #632 

FO  / 6092. 84/,  XK  /273.15/ 

BIRD  CALIBRATION  DATA  FOR  SENSOR  #638  —  9-18-87  CAL. 


SEA 

DATA 

DATA 

SEA 


DATA  A  /3.67399029E-03/ ,  B  /6.01229226E-04/ ,  !  SEABIRD  SENSOR 

1  C  /I . 51578463E-05/ ,  D  /2.71329547E-06/  !  #638 

DATA  FO  /6337.28/,  XK  /273.15/ 


f 


I 


CALL  HEXDEC(FRQ(1)) 


C 

C*** 

c 


IF(FREQ.NE.O. )  THEN 
DIV=ALOG(FO/FREQ) 
T1=A+B*DIV 
T2=C*DIV**2 
T3=D*DIV**3 
TEMP=1./(T1+T2+T3)-XK 

ELSE 

TEMP=0. 

END  IF 

END  OF  SUBROUTINE  CNVTMP 

RETURN 

END 


:i  ^5 
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c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  CNVCND 

CHARACTER*6  FRQ(8) 

CHARACTERS  ANLG(8) 

REAL  IEXP 

COMMON  /  RAW  /  ANLG , FRQ 

COMMON  /  CNVRT  /  C0ND,SAL78, TEMP, VOLT, FREQ, XL4,XL8,PSI, 

1  DEPTH, TRANS, TLIGHT,CHLA 

SEA  BIRD  CALIBRATION  DATA  SENSOR  #267 

DATA  W  /4.43054668E-10/,  X  /5.20135966E-01/,  !  SEA  BIRD  SENSOR 

1  Y  /-4. 36199682/,  Z  /2.42254149E-04/  t  SENSOR  #267 

DATA  IEXP  /8.4 / 

SEA  BIRD  CALIBRATION  DATA  SENSOR  #234  —  9-18-87  CAL. 


1 


DATA  W  /8.88229858E-09/,  X  /5.1687217E-01/, 
Y  /-4. 37109279/,  Z  /-4.80212277E-05/ 
DATA  IEXP  /7.2/ 


SEA  BIRD  SENSOR 
!  SENSOR 


#234 


CALL  HEXDEC(FRQ(2)) 


F=FREQ/1000. 

C1=W*F**IEXP 

C2=X*F**2 

C3=Y 

C4=Z*TEMP 

C0ND=Cl+C2+C3+C4 


P 


CALL  SALNTY 


l 


P  ** 


END  OF  SUBROUTINE  CNVCND 

RETURN 

END 


> 


I 


SUBROUTINE  SALNTY 

COMMON  /  CNVRT  /  COND,SAL78, TEMP, VOLT, FREQ, XL4,XL8,PSI, 
DEPTH , TRANS , TLIGHT , CH  LA 

T.TEMP 

P=*2.  !  HARD  CODED  FISH  DEPTH  FOR  SAL  CALC. 

SAL78=0. 


>  IF(COND.GE. 0.0005)  THEN 

DT=T-15. 

R=C0ND/42.914 

RT35=((( . 10031E-08*T-.69698E-06)*T+. 0001 104259 )*T 
1  +. 200564E-01)  *  T+. 6766097 

RT=R/(RT35*(1 . +PNC(P)/(FNB(T)+FNA(T)*R) ) ) 

I  RT=SQRT(ABS(RT) ) 
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XR=RT 


1 


XT=DT 

Sl=( ( ( ( 2 . 7081*XR-7 . 0261 )*XR+14 . 0941 )*XR+25 . 3851 )*XR- . 1692 )*XR 
S2=.008 

S3=XT/(1.+.0162*XT) 

S4= ((((-. 01 44*XR+ . 0636 ) *XR- . 037 5 ) *XR- . 0066 ) *XR- . 0056 ) 
*XR+.0005 
SAL78=S1+S2+S3*S4 
END  IF 


C 

C*** 

C 


END  OF  SUBROUTINE  SALNTY 


RETURN 

END 


SUBROUTINE  CNVPSI 


CHARACTER*4  ANLG(8) 

CHARACTER*6  FRQ(8) 

COMMON  /  RAW  /  ANLG.FRQ 

COMMON  /  CNVRT  /  COND  f  SAL78 ,  TEMP ,  VOLT ,  FREQ  ,"1.4,  XL8 ,  PSI , 
1  DEPTH, TRANS, TLIGHT,CHLA 


1 


DATA  C  /-57083.99/,  D  /O. 03636303/,  TO  /27. 597905/, 
G  /6.894759E7/,  H  /10.0E-6/,  XLAT  /36.0/ 


CALL  HEXDEC(FRQ(3)) 

F=FREQ 

PSI=G*(C*( 1 . -(H*T0*F)**2) )*( 1 . -D*( 1 . -(H*T0*F)**2) ) 


X=SIN(XLAT/57. 29578) 
X=X**2 


GRAVTY=9 . 780318*( 1 . + (5 . 2788E-3+2 . 36E-5*X)*X)+1 . 092E-6*PSI 
DEPTH=( ( (-1 . 82E-15*PSI+2 . 279E-10)*PSI-2 . 2512E-5)*PSI+9 . 72659)*PS 
DEPTH=DEPTH/GRAVTY 


I 


RETURN 

END 


C 

C 

C*** 

C 


FUNCTION  FNC(XP) 

FNC=( ( . 3989E-14*XP- . 637E-09)*XP+ . 207E-04)*XP 

END  OF  FUNCTION  FNC 

RETURN 

END 
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c 

c*** 

c 


- *  —  —  —  — # - •  —  —  —  —  • 

FUNCTION  FNB(XT) 

FNB= ( . 4464E-03*XT+ . 03426 ) *XT+ 1 . 

END  OF  FUNCTION  FNB 

RETURN 

END 


F 

c 

c*** 

c 


FUNCTION  FNA(XT) 

FNA=- . 003107*XT+ .4215 

END  OF  FUNCTION  FNA 

RETURN 

END 


SUBROUTINE  CNWLT(AN) 


CHARACTERS  0,AN 
CHARACTER* 1  NUM 
INTEGER*2  G 

COMMON  /  CNVRT  /  COND.SAL78, TEMP, VOLT, FREQ, XL4,XL8,PSI, 
1  DEPTH, TRANS, TLIGHT,CHLA 


P 


REFERENCE  VOLTAGE  —  VREF  —  FOR  TUPS  COMPUTER 


P 


DATA  VREF  71.23536/ 


ONE  COMPUTER  OLD 


DATA  VREF  /1.233/  !  OTHER  COMPUTER  NEW  —  USED  6/88 

0  =  AN 

LENGTHS 

XN=0. 

P-1. 

!  DO  100  I=LENGTH, 2 , -1 
j  G  =  ICHAR(0(I : I ) ) 

!  IF(G.GT.64)  THEN 

i  D=FLOAT(G)-55 . 

ELSE 

NUM  =  AN(I : I ) 

IF  (NUM.EQ. '0' )  D  =  0.  !  ADDED  BY  ARNONE  2-5-1988 

IF  (NUM.EQ. '1')  D  =  1. 

IF  (NUM.EQ. '2')  D  =  2. 

IF  (NUM.EQ. '3')  D  =  3. 

IF  (NUM.EQ. '4')  D  =  4. 

IF  (NUM.EQ. '5')  D  =  5. 

IF  (NUM.EQ. '6')  D  =  6. 

IF  (NUM.EQ. '7')  D  =  7. 

IF  (NUM.EQ. '8')  D  =  8. 

IF  (NUM.EQ. '9')  D  =  9. 

END  IF 


XN=XN+D*P 
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_  •  _  —  _ • _ _ • •  _____  __  •  _  _ _  _  •  

- • •  —  • •  —  —  • - • •  •  •  •  •  •  • 

P=P*16. 

100  CONTINUE 

V0LT=XN*VREF/2048 . 

IF  (AN(1 : 1) .EQ. 'C' )  V0LT=-V0LT 

o* 

mo 

:***  END  OF  SUBROUTINE  CNWLT 

RETURN 

END 


SUBROUTINE  DIVIDE(RRAW) 


SUBROUTINE  DIVIDE  SEPARATES  THE  INPUT  STRING  FROM  TUPS 
IN  10  FREQUENCY  AND  ANALOG  CHANNEL  OUTPUT 

VARIABLES: 

FRQ  -  ARRAY  TO  HOLD  FREQUENCY  DEFINING  STRINGS 
AN LG  -  ARRAY  TO  HOLD  ANALOG  DEFINING  STRINGS 


CALLING  PROGRAM(S)  : 
PROGRAM  CNVRT 


CHARACTERS  *)  RRAW 
CHARACTER**  ANLG(8) 
CHARACTER*6  FRQ(8) 
CHARACTER*80  0UTSTR(12) 
LOGICAL* 1  I LOG (1600) 

EQUIVALENCE  ( ILOG( 1 ) f  OUTSTR) 

INTEGER  LSTR(20) 


COMMON  /  RAW  /  AN LG , FRQ 

CALL  PARSE (RRAV, OUTSTR, NSTR) 

NFRQ  =  0 
NANAL=  0 
DO  I  =  1 ,NSTR 

LSTR(I)  =  INDEX (OUTSTR (I),'  ')-l 
IF(LSTR(I)  .GT.  1)  THEN 
IF(I.GT. 1  .AND.  LSTR(I-l)  .EQ.  1  )THEN 
NFRQ  =  NFRQ  +  1 

FRQ(NFRQ)  =  0UTSTR(I-1)(1: 1)//'  ' //OUTSTR(I ) ( 1 : 4) 
ELSE 

NANAL  =  NANAL  +  1 
ANLG(NANAL)  =  0UTSTR(I)(1:4) 

END  IF 
END  IF 

END  DO 
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RETURN 

END 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


SUBM 


OUTINE 

IPOSE 


PUR 


PARAMETERS 


PARSE( INSTRING , OUTSTR , N) 

BREAKS  A  STRING  INTO  AN  ARRAY  OF  CHARACTER  STRINGS 
WITH  EACH  BEING  ONE  ELEMENT  OF  A  LIST.  LIST 
SEPARATORS  ARE  COMMAS  AND/OR  BLANKS 

INSTRING  -  A  CHARACTER  STRING 

OUTSTR  -  AN  ARRAY  OF  CHARACTER*80  BY  10 

N  -  NUMBER  OF  SUBSTRINGS  FOUND 


SUBROUTINE  PARSE( INSTRING, OUTSTR, N) 


CHARACTER* ( * )  INSTRING, OUTSTR(*) 

INTEGER  N , PTR , DPTR 

LOGICAL  MORE , BLANK , COMMA 


MORE  =  .TRUE. 
PTR  =  0 

N  =  0 

DO  WHILE (  MORE  ) 


IBLANK  =  INDEX( INSTRING (PTR+1:),'  ') 
ICOMMA  =  INDEX (INSTRING (PTR+1 
I  QUOTE  =  INDEX(INSTRING(PTR+1 

DPTR  =  MIN(IBLANK, ICOMMA) 

IF( IBLANK. EQ.O)  DPTR  =ICOMMA 
IF ( ICOMMA. EQ.O)  DPTR  =IBLANK 


C 

C 


TYPE  * , ' IQUOTE  =  ', IQUOTE 

TYPE  * , ' DPTR , IBLANK , ICOMMA= ' , DPTR , IBLANK , ICOMMA 


IF( IQUOTE  .NE.  0  .AND.  IQUOTE  .LT.  DPTR  )THEN 
NQU0TE=INDEX(INSTRING(PTR+IQU0TE+1 s ) , ' '  "  )+IQUOTE 


E 


E 


IF(NQUOTE.EQ. IQUOTE)  THEN 
IQU0TE=0 
ELSE 

IBLANK=INDEX(INSTRING(PTR+NQU0TE+1: ) , '  ' )+NQUOTE 
ICOMMA=INDEX( INSTRING (PTR+NQUOTE+1 
IF( ICOMMA  .NE.  0)  IC0MMA=IC0MMA4NQU0TE 
END  IF 

TYPE  * , ' IQUOTE , NQUOTE , IBLANK , ICOMMA , PTR ' 

TYPE  * , IQUOTE , NQUOTE , IBLANK , ICOMMA , PTR 
ELSE 

IQUOTE  =  0 

TYPE  *, 'I QUOTE, NQUOTE, I BLANK, I COMMA, PTR' 

TYPE  * , IQUOTE , NQUOTE , IBLANK , ICOMMA , PTR 
END  IF 


BLANK  =  .TRUE. 

DO  I  =  PTR+1, 80 

IF(INSTRING(I : I )  .NE.  '  ')  THEN 
BLANK=  .FALSE. 
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GOTO  1000 

END  IF 

END  DO 
CONTINUE 

IF(BLANK)  IBLANK=0 

IF(IBLANK  .EQ.  0  .AND.  ICOHMA  .EQ.  0)  THEN 
MORE  =  .FALSE, 

ELSE 

COMMA  =  .FALSE. 

DPTR  =  MIN(IBLANK,ICOMMA) 

IF(IBLANK.EQ.O)  DPTR  -ICOMMA 
IF(ICOMMA.EQ.O)  DPTR  =IBLANK 
IF(ICOMMA  .EQ.  DPTR)  COMMA  =  .TRUE. 

N  =  N  +  1 

IF(COMMA  .AND.  DPTR  .EQ.  1)  THEN 
OUTSTR(N)  *  '  ' 

PTR=PTR+1 

ELSE 

OUTSTR(N)=INSTRING( PTR+1 : PTR+DPTR-1 ) 


IF(IQUOTE.NE.O)  THEN 

OUTSTR(N) =INSTRING ( PTR+l+IQUOTE : PTR+NQUOTE-1 ) 
END  IF 


P 


TYPE  * , N, I BLANK , ICOMMA , OUTSTR(N) 
PTR=PTR+DPTR 

IF(OUTSTR(N)  .EQ.  '  ')  N=N-1 
END  IF 

END  IF 

END  DO 
RETURN 
END 


30 


subroutine  time£ix(time, timea,n) 
real*8  time(n) , timea(n) 
real*8  dt 
inum=n 
do  i=l,15 

type  *,i,time(i) 
if(time(i).eq. time(i+l)  )  then 
else 

kks=i+l 
goto  30 

endif 

enddo 
continue 
do  ii=l , 15 

i=inum-ii+l 
type  *,i,time(i) 
if(time(i).eq. time(i-l))  then 
else 

kke=i 
goto  31 

endif 

enddo 

nt  =  kke  -  kks 

dt=  (time(kke)  -  time(kks))/float(nt) 
type  *,nt,dt 
timea(kks)=time(kks) 
do  i=kks-l,l,-l 

timea( i )=timea( i+l)-dt 
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enddo 

do  i=kks+l, inum 

timea(i)=timea(i-l)+dt 

enddo 

return 

end 


sub rou tine  sec  to  hms (ts,ih,im,is) 
implicit  real*5  (t) 
ih=ts/3600. 
tl=ts-ih*3600. 
im=tl/60. 
tl=tl-im*60. 
is=tl 

return 

end 


UOOOUUOOUUOOOUO 
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PR0|G)RAM  READJUN 

PUR|P|0SE  1)  program  readjun  to  deal  with  tess  files 

2)  place  tess  file  in  appropiate  structure 

AtnUoR  NORMAN  GUINASSO 

DAT®  JUNE  88 

DEC.  88 


REVffl 


SION 


implicit  real*8  (d) 

real*8  ttime,  flatd,  flatm,  flats ,  flond, flom, flons 
real*8  d60  /60.d0/ 

REAL*4  V(ll),spd 

CHARACTER*21  DATELINE 
CHARACTER*80  dataline 
character*80  naviline 
character*40  logname 
character*20  cdate 

type  '("  Program  to  read  and  translate  TESS  files  ")' 
type  '("  Enter  input  file  name  =>  ",$) 
accept  ' (a)' , logname 

cdate='  '  !  put  date  and  time  into  string  cdate 

call  date(cdate) 
call  time(cdate(12: )) 

!  ope(n|  files 

open ( 9 , name=logname , readonly , s  t a  t us= ' old ' , err«9^9 ) 
call  chg_ext( logname, logname, ' rdj ' ) 

open( 11 , carriagecont rol= 'list' , s tatus=' new' , name=logname , recl=255|) 
type  Output  file  will  be  called  ' ' ,a) ' , logname 


write(ll, '  ("TESS  data  processed  by  program  r  djun  on  ",a)')cda)t 
vrite(ll, ' (' 'TESS  data  are  2.5VREF,  441  NM,  4b8  NH,  PYR01,PYR02 
FILENAME  =  "  ,a)' ) logname 
icnt=0 

do  while  ( . true. ) 

read(9, ' (a)' , iostat=ios)  dateline 
if(ios.eq.-l)  goto  100  !  if  end  of  file 

icnt=icnt+l 

j=str$translate(dateline, dateline, ' 
read (dateline,*, iostat=jos)  mon, iday , iyear , ih, im, is 
if(jos.ne.O)  type  *,'  read  error  on  dateline  string  ',jos| 
if  (iyear. gt .  1900)  iyear=iyear-1900  !  limit  to  20th  centujry 
ttime=d60*d60*(ih  +(is/d60  +im)/d60)!  seconds  since  midni 


%ht 

if (mod(icnt ,50) ,eq .  1)  type  '("  +  ", i5, 2X) ', icnt !  preveNt  Iboredom 

read(9, ' (a) ' , iostat=ios)  dataline 
read(9,'(a)' ,iostat=ios)  naviline 

j=str$translate(naviline, naviline, '  ','NV')  !  eliminate  tlW 

read(naviline,*, ios tat=jos) flatd , flatm, flats, flond, flonm, flons, 
ihdg, spd, tdl , td2  !  get  the  numbers 

if(jos.eq.-l)  goto  110  !  if  they  are  not  there,  stop 

if(jos.ne.O)  type  *,'  error  reading  naviline  ',jos!  signal 
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|100 

110 

115 

120 

999 

300 


dlat=  flatd  +(flats/d60  +  flatm)/d60 
dlon=  flond  +(flons/d60  +  flonm)/d60 
dlon=  -dlon 

j =str$trim(dataline, dataline, len) 


!  decimal  degrees) 
!  ditto 
!  Occident 
I  remove  trailing) 


ROUTINE  TO  CONVERT  VOLTS  IN  THE  DATA  LINE  TO  ENGINEERING  UNITS 
READ( DATALINE , * , IOSTAT=KOS )<V(K),K=1,11) 
if (ROS.eq.-l)  goto  115  !  if  they  are  not  there,  stop 

CALL  CALCENGR(V) 

vrite(ll,300)  iyear,mon, iday , ttime, ih, im, is,  !  record  for) 

dlat,dlon,tdl,td2,ihdg,spd,V(2),(V(K),K=5,8)  !  posterity 

enddo 

type  *,icnt,  '  records' 
stop  'end  of  file' 

stop  'unexpected  end  of  file  reading  dataline' 

stop  'unexpected  end  of  file  reading  dataline  FOR  CALCULATION' 

stop  'unexpected  end  of  tile  reading  naviline' 

stop  'bad  inputfile' 

format(i2.2, '/' , i2.2, '/' , i2.2, lx, f8. l,2x, i2.2, ' j ' , i2.2, ' : ' , i2.2, l)x, 
fll.6,fl2.6,2f9.2,2x,i3.3,F5.1,5F7.3) 

end 


blanks 


i 


SUBROUTINE  CALCENGR(V) 

REAL*4  V(ll) 

SUBROUTINE  TO  CONVERT  TESS  VOLTAGES  TO  ENGINEERING  UNITS 
THE  VOLTAGE  ARRAY  "V"  IS  PASSED  BACK  AND  FORTH 


> 


C 

C 


V(5)  IS  441  LOG  LIGHT  SENSOR  SNs  7102 
UNITS  ARE  uW/cm2  sec  nm 
V(5)=3.54*EXP(V(5)-0.734) 

V(6)  IS  488  LINEAR  LIGHT  SENSOR  SNs  7105 
UNITS  ARE  uW/cm2  sec  nm 
V(6)=V(6)/0. 25599 


P 

P 


> 


V(7)  IS  PYR01  —  THE  LIGHT  BULB  PYROHELI 0METER 
UNITS  ARE  uV/cm2*10"4 

V(7)=(V(7)-0.0)/0.640  !  CONVERT  VOLTS  TO  MILLIVOLTS 

V(7)-V(7)*(0. 1202/60. 0)*4. 184*100.0  !  CALIBRATION  AND  CONVERSION 


► 


C 

C 


V(8)  IS  PYR02  —  THE  HEMISPHERE  PYROHELI OMETER 
UNITS  ARE  uW/cm2*10“4 

V(8)=(V(8)-0.0)/0.9719  !  CONVERT  VOLTS  TO  MILLIVOLTS 

V(8)=V(8)*(0. 2232/60. 0)*4. 184*100.0  !  CALIBRATION  AND  CONVERSION 


fc** 


END  OF  SUBROUTINE  CALCENGR 
RETURN 
END 


i 


> 


> 


» 


I 


uoooooouuoooou 
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|c*** 

PR 

PURPOSE 


.OjGJRAI 


AUTtHOR 
DA' 
REVtUSION 


*************************************************************** 
M  PRETABLE 

1)  reads  in  tups  "cnv"  file 

2)  filters  some  data  and  writes  a  binary  file 
NORMAN  GUINASSO 
JUNE  88 
DEC.  88 


Q****)*)***************:*****:*  Air  ***★*★***★★*★★★★*:★******★*★★*****•**★*★★****+* 

c 


! 


da 


parameter  (nn=5000) 
common  /fopen  name/  f  name 
charac t er*60  T_name , filename , ou t  f i le 
tatups  variables: 

equivalence  (datatups(l,  1),  tetnp(l)) 
equivalence  (datatups(l,  2),  salin(l)) 
equivalence  (datatups(l,  3),  ptrans(l)) 
equivalence  (datatups(l,  4),  fluor(l)) 
equivalence  (datatups(l,  5),  d465nm(l)) 
equivalence  (datatups(l,  6),  d507nm(l>) 
equivalence  (datatups(l,  7),  d532nm(l)) 
equivalence  (datatups(l,  8),  transv(l)) 
equivalence  (datatups(l,  9),  flsigv(l)) 
equivalence  (datatups(l,10),  flsclv(l)) 
equivalence  (datatups^lfll),  d488up(l)) 
equivalence  (datatups(lf 12),  dep(l)) 


real*8  datatups(nn, 12) 

real*8  temp(nn) , salin(nn) , ptrans(nn) , transv(nn) 
real*8  fluor(nn) , flsigv(nn) , flsclv(nn) 
real*8  d465nm(nn) ,d507nm(nn) fd532nm(nn) ,d488up(nn) 
real*8  dep(nn) 


real*8  ravdepth(nn) 
real*8  rawfluor(nn) 
real*8  ddint(ll) 
real*8  data(nn,ll) 
real*8  ttime(nn) 
real*8  out(nn) 
character*150  in 
character*8  ctime, timestr 
character  ans 

character*30  timedate, f iledate 


call  date(timedate(ll: )) 
call  time( timedate(21: )) 
timedate(l : 10)=  'Processed' 


type  *,  '  Preprocessor  for  TUPS  file  ' 

type  *,  '  This  program  filters  fluorescence,  X  trans,  depth,  sal 
type  *  ,'  and  writes  a  binary  file' 
call  fopenread(lO) 
type  *,  '  ' 


inity' 


i=0 


do  while( . true. ) 

read(10, ' (a) ' , iostat=ios)  in 
if(ios.eq.-l)  goto  290 


CL  CL  Q.  Q. 
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i=i+l 

if (mod(i ,  100) .eq.  1)  type  '("  +  .",$)' 

j=str$translate(in,in, '  ',':/|') 

j=str$trim  (in,in,len) 

read ( in( 1 : len )  ,*, ios  tat =ios ) 

iyr, imo, iday, ttime(i) , ihr, imin, isec, 

(datatups(i ,j ) , j=l, 12) 
if(i.gt.l)  then 

if(  ttime(i)  .It.  ttime(i-l)  )  then 
i = i  —  1 

type  *,'  bad  time  at  i,ttime(i-l) 
type  *,in(l:70) 

endif 

endif 


C90 


! 

104 


! 


t 

i 


if(ios.ne.O)  then 
i=i-l 

type*, ios,'  ',in(l:66) 

else 

endif 

enddo 
continue 
ncnt=i 
do  i=l,ncnt 

rawdep t h ( i ) =dep ( i ) 
ravfluor(i)=fluor(i) 

enddo 

write(filedate,104)  iyr, imo, iday 
format('filedate  ' ,i2, ' /' ,i2.2, ' /' ,12.2) 


type  *,ncnt 

type  *,  '  salinity  ' 


call  rfilter(salin,ncnt,  .03,  2.1  , 

pause 

30., 

37.) 

type  *,  '  percent  trans' 

call  rfilter(ptrans,ncnt,  .1,  6.  , 

60., 

iOO.) 

pause 

type  *,  '  depth  ' 

call  rfilter(dep,ncnt,  .3,  5.  , 

o., 

30.) 

pause 


type  *,  'fluorescence' 
f=.5 

fmin=100. 
f max-850. 

call  rfilter  (  fluor,  ncnt,  f,  65.,  fmin,  fmax  ) 


open(15,name='plotfluor.out' ,status='nev' ,carriagecontrol='list' ) 
write(15,'(f7.3,3f7.1)') 

#  ( ttime(i)/3600. , ravfluor(i) ,rawdepth(i) , fluor(i) , i=l,ncnt) 

close(15) 

call  chg_ext(f_name, out file,  'pto'  ) 


open(unit=14,name=outf ile,status='new' , iostat=jos, 
form=' unformat  ted' ,carriagecontrol='none' ) 
type  *, jos, '  '//outfile 
write  (14)  timedate, f iledate,ncnt 
type  *,  '  ncnt  ',ncnt 
do  i=l,ncnt 

vrite(14)  t timer' ) , (datatups(i , j ) , j=l, 12) 

enddo 

close(14) 
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!  wine 
!  if  r 
!  als( 
!  guii 


stop 

format( '  +  . ' ,$) 

format(6i3.2,lx,a,lx,f8.1,3(/,8£10.2)) 

format(f9.4,2f8. 1) 

end 

subroutine  rf ilter(d,n, f ,windovmin, fmin, fmax) 
swing  filter 

sw  point  differs  by  previous  point  times  f  then  old  point  is  retained 
if  point  out  of  range  fmin, fmax  than  previous  point  is  retained 
isso,  circa  1988 
real*8  f,vindovmin, fmin, fmax 
real*8  d(n) 
real*4  ave  ,  sdev 

type  *,  '  rfilter-i-f=  ' , f ,vindovmin 

type  *,  '  rfilter-i-fmin,fmax=  ', fmin, fmax 

call  stat(d,n, ave, sdev) 

i changed  =  0 

tml=0 

tm2=0 

do  i=2,n 

tm3=tm2 

tm2=tml 

tml=t 

t=d(i) 

prev  =  d(i-l) 

if(prev.lt.fmin  .or.  prev. gt. fmax)  then 
prev  =  ave 

type  prev  replaced  with  average  =  ' ,ave 

end  if 

delta  =  abs(d(i)  -  prev) 
window  =  f  *  prev 

if (window. It. windowmin)  window=windowmin 

if (delta. gt. window  .or.  d(i).lt.fmin  .or.  d(i).gt.fmax)  then 
d(i)  =  prev 
TYPE  ' (Ix,2f8.2,4f8. 1, 

"  changed  ",i5,"  to  ",f8.1)', 
window, del ta,  tm3,  tm2,’  tml ,  t ,  i ,d(i) 
i changed  =  i changed  +  1 

else 
end  if 

enddo 

return 

end 


subroutine  table_lookup  (x,y,n,nn,m,xx,yy, reset) 

linear  interpolation  lookup  program 
guinasso  circa  1988 

n  dimension  of  x,y 

nn  number  actually  used 


real*8  x(n) 
real*8  y(n,m) 
real*8  xx 
real*8  yy(m) 
real*8  f 

integer*^  il  /l/ 
logical  reset 


array  of  x's  must  be  monotonically  increasing 
array  of  y's 

x  value  at  point  to  look  up 

y's  at  xx  (returned) 

interpolation  fraction 

start  of  interpolation  range 

will  reset  range  looked  at  to  1  thru  n 


|  if  (reset)  i 1= 1 
j  if(xx.lt.  x(l)  )  goto  10 
do  i=il,n-l 

j!  if(x(i).ge.x(i+l))  then 


!  not  on  table 
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stop  'table  lookup-f-bad  table,  not  monotonic' 
else  if(xx.ge.x(i)  .and.xx.le.x(i+l))  then 
nl=i 
n2=i+l 
il=nl 
goto  20 

end  if 


enddo 

continue 


type  *,  'table  lookup-i-xx  not  on  table  ' 
type  *,  xx,x(T),x(n) 


return 


£  =  (xx-x(nl)  )/(  x(n2)-x(nl)  ) 
wf  =  l.dO-f 
do  j=l,m 

yy(j)  =  vf*y(nl, j )  +  f*y(n2,j) 

enddo 

return 

end 

character*8  function  timestr(secs) 


!  guir 
!  -  - 


»rts  seconds  past  midnight  to  HH:MM:SS 
isso  circa  1988 


ih=secs/3600 

im=(secs  -  ih  *  3600. )/60. 
is=nint(secs  -ih  *  3600.  -  im  *  60.) 

write  (timestr, ' (i2.2, "  :  "  ,12.2, "  : "  ,i2.2)'  )  ih,im,is 

return 

end 

subroutine  stat(y,n,aver,sdev) 
llates  mean  and  standard  deviation  of  array  y 
real*8  y(n) 
sum=0. 

do  i=l,n 

sum=sum+y(i) 

enddo 

aver  =  sum/float(n) 

sum=0. 

do  i=l,n 

sum=sum+  (aver-y(i))**2 

enddo 

var  =  sum/float(n-l) 
sdev  =  sqrt(var) 

type  *, 'stat-i-mean,  standard  deviation  =  ',  aver, sdev 
return 


! 


end 
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>*** ***************************************************************** 

«/ 

:  PROGRAM  MERGETT 

-/ 

:  PURPOSE  1)  merge  norda  tess  tups  files 

2  AUTHOR  NORMAN  GUINASSO 

L/ 

G  DATE  JUNE  88 

r% 

:  DEC.  88 


SUBROUTINES 


****************************************************************** 

smical  and  Environmental  Resaearch  Group 
s  A&M  University 

mxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

parameter  (nn=5000) 
parameter  (nnn=25000) 

common  /fopen  name/  f_name 
character*60  T  name,outf ile, infile 
character*150  In 
real*8  ddint(ll) 
real*8  data(nn,ll) 
real*8  datatups(nnn, 12) 
real*8  attime(nnn) 
real*8  time(nn) , ttime 
character*8  ctime, timestr 
character*30  timedate, f iledate 

call  lib$init_timer() 
type  *,  '  tess  file  ' 
call  fopenread(9) 


tups  file 

call  chg_ext(f  name, inf ile, ' pto' ) 

type  Opening  TUPS  pretable  output  file  infile 

open  (14,  name=inf ile, form='unformatted' ,status='old' ,iostat= 

if(jos.ne.O)  stop  'file  not  found' 

put  file 

call  chg_ext(f_name, out file, ' fin' ) 
type  *,'  Output  file  will  be  called  ',outfile 
type  *, '  ' 
tess  file 

do  vhile( . true. ) 

read(9, ' (a) ' , iostat=ios)  in 
if (ios.eq.-l)  goto  90 
i  =  i  +  l 

j=str$translate(in, in, ' 
read(in,*, iostat=ios) 

*  iyr, imo, iday, time(i) , ih, im, is, (data(i , j ) , j=l, 11) 

if(ios.ne.O)  then 

type  *,in(l:76) 
i=i-l 

end  if 

i  if (mod( i , 100) .eq . 1)  type  100 


enddo 
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JO  nctess=i 
close(9) 

type  *,  nctess,  '  records  in  tess  file  ' 
type  *,time(l),time(nctess),  '  are  start  stop  times' 
call  lib$shov_timer  (  ) 
read  (14)  timedate, filedate, nctup 
type  *,  nctup,'  ' , timedate, filedate 
do  i=l, nctup 
k  =  i 

read(14,end=92)  attime(i) , (datatups(i , j ), j=l,12) 

enddo 
goto  93 

J2  type  end  of  tups  file  ',k, nctup 

if (k.ne. nctup)  type  *,  'bad  tups  file' 

J3  close(14) 

call  lib$show  timer() 


type  starting  output  phase' 

type  ' 

open(15,name=outfile,carriagecontrol=' list' ,status='new' , 

*  recl=256) 
iout=0 

do  1=1, nctup 

if(mod(i,100).eq.0)  type  100 
ttime  =  attime(i) 

call  table  lookup(time, data, nn, nctess, 11, 

*  t time, ddint, . false. ) 
ctime=timestr( ttime) 
iout=iout+l 

vrite(15,101)iyr,imo,iday,ctime, 

*  ttime, ddint, (datatups(i,j),j=l,12) 
enddo 

call  lib$show_timer() 
stop 

100  format('+. ' ,$) 

101  format(3i3.2,lx,a,lx,f8.1,fl0.6,fll.6,2fl0.2,f6.1,f5.2, 

*  5f7 .3,  f6.2,  f7.3,  f6.2,  f9.2,  7f7.3,  f6.2) 

!  tess 

!  F  la t, long, tdl , td2, heading, speed, 2. 5  ref, data  1-4 
!  S  temp,  sal,  ptrans,  chlorfl,  data  1-7,  depth 
end 

! 

subroutine  table  lookup  (x,y,n,nn,m,xx,yy, reset) 

! - - - 

!  linear  interpolation  lookup  program 
!  guinasso  circa  1988 


real*8  x(n) 
real*8  y(n,m) 
real*8  xx 
real*8  yy(m) 
real*8  f 

integer*4  il  111 
logical  reset 


dimension  of  x,y 
number  actually  used 


array  of  x's 
array  of  y's 
x  to  look  up 
y's  at  xx 

interpolation  fraction 
start  of  interpolation  range 
will  reset  range  looked  at  to 


if  (reset)  i 1=1 

if(xx.lt.  x(l)  .or.  xx.gt.x(nn))  goto  10 
do  i=il,nn-l 

if(x(i).ge.x(i+l))  then 

type  *,  i,x(i),x(i+l) 

stop  'table  lookup-f-bad  table' 


thru  n 


!  not  on  table 
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end  if 

if(xx.ge.x(i)  .and.xx.l=.x(i+l))  then 
nl=i 
n2=i+l 
il=nl 


10 


goto  20 

end  if 

enddo 

continue 

type  *,  ' table_lookup-i-xx  not  on  table  ' 
type  *,  xx,x(l),x(nn) 

return 


! 

20 


f  =  (xx-x(nl)  )/(  x(n2)-x(nl)  ) 
wf  =  l.dO-f 
do  j=l, m 

yy(j)  =  wf*y(nl, j )  +  f*y(n2,j) 

enddo 

return 

end 


!  - 


!  converts  seconds  past  midnight  to  HH:MM:SS 
!  guin|asso  circa  1988 
!  -  - 


character*8  function  timestr(secs) 


ih=secs/3600 
im=(secs-ih*3600. )/60. 
is=nint(secs-ih*3600.-im*60. ) 

write  (timestr,'(i2.2,":",i2.2,":",i2.2)'  )  ih,im,is 

return 

end 


ooooouuuouou 
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c*** 

PROjGjRA! 
PUR 
Ainwo 
DA' 

REVl 


vm 


*************************************************************** 
M  TTPLOT 

1)  plots  combined  tess  and  tups  data 
NORMAN  GUINASSO 
JUNE  88 
DEC.  88 


|P(0SE 
R 


SION 


|q*  *★★★★★★*★★★★***★★**★*★  ’k'k'k'k'k'k'kic'k'kirkirk'k'k'k’kirk'kirk'k’k'k'k'kirk'k'k’kicic'k'k'k'kitirki^k 

parameter(nn=5000) 
common  /paper/  xpaper,ypaper, LASER 
common  /xrange/  xl,x2 
LOGICAL  LASER  /.FALSE./ 
real*4  time(nn) 
real*4  flat(nn) 
real*4  flon(nn) 
real*4  hdg(nn) 
real*4  speed(nn) 

real*4  fli tel(nn) , fli te2(nn) , flite3(nn) , £lite4(nn) 
real*4  temp(nn) , salin(nn) 
real*4  ptrans(nn) 
real*4  fluor(nn) 
real*4  depth(nn) 

real*4  upl(nn) ,up2(nn) ,up3(nn) ,downl(nn) 
character*4  itime 
CHARACTER*300  IN 

character*60  tl, t2, t3, t4, t5, t6, t7, t8 
character*30  pfile 
character*30  filename 
character*2  cday 
data  z,sx  /0. ,42./ 

______  -1234567890  - - - 

f ilename='xxjun88. f in' 

tl=' Geochemical  and  Environmental  Research  Group' 
t2='Texas  A&M  University' 
t3='Ten  South  Graham  Road' 
t4='College  Station,  Texas  77840' 
t5='Telephone:  409  690  0095' 
t6='N0RDA  Continuous  Underway  Data' 
t7='  ' 

t8=' Department  of  Oceanography' 
j=str$trim(t6, t6,len) 
type  '("  enter  day  of  plot  ##> 
accept  '(a)', cday 
t7=cday//'  June  1988' 
t6=t6( 1 : len+l)//t7 
type  *, t6 
type  *,t7 
k=0 

£ilename(l:2)=  cday 

open(9,name=filename,status='old' , iostat=ios) 
type  *,  'iost  t  =',ios,'  '//filename 
if(ios.ne.O)  stop 

DO  VHILE( .TRUE. ) 

READ(9 , ' (A) ' , END=50)  IN 
k=k+l 

j=str$trim(in,in,ilen) 
if (mod(k, 100) .eq. 1)  type 
j=str$translate( in, in, ' 


("  +  . 


,$)' 
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t 
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if(k.lt.lO)  then 

type  *,'$l'//in(l:150) 

if (ilen.gt . 150)  type  *, '$2'//in(151 : ilen) 

end  if 

read(in,*,iostat=ios)  iyr,imo, iday, ihr, imin, isec, 
time(k) , flat(k) , f lon(k) ,  tdl , td2 ,hdg(k) , speed (k) , ref , 
f li tel(k) , f lite2(k) ,flite3(k), flite4(k) , temp(k) ,salin(k) , 
ptrans(k) , fluor(k) ,upl(k) ,up2(k) ,up3(k),d4,d5,d7,downl(k) ,depth(k) 

type  *, iyr, imo, iday , ihr , imin, isec, 
time(k) , flat(k) , f lon(k) , tdl , td2,hdg(k) , speed (k) , ref, 
flitel(k),flite2(k),flite3(k),flite4(k), temp(k) ,salin(k) , 
ptrans(k) , fluor(k) ,upl(k) ,up2(k) ,up3(k) ,d4,d5, d7, down! (k) , depth (k) 

if(k.gt.250)  goto  60 

enddo 

continue 

continue 

pfile  =  ' ttplotxxZ. iop' 
pf ile(7:8)=cday 
IF(LASER)  THEN 

pfile(9:9)='L' 

ELSE 

pfile(9:9)='P' 

ENDIF 

type  *, 'starting  plot  '//pfile 
xpaper=41.50 
ypaper=32.5 
call  iopen(20, pfile) 

IF(  LASER) call  set_42  (z,sx,z,sx,  z,sx,z,sx,l) 

IF( . NOT . LASER) call  set_60  (z,sx,z,sx,  z,sx,z,sx,l) 
n=k 

do  i=lfn 

time(i)=time(i)/3600. 

enddo 

ixl  =  time(l) 
ix2  =  time(n)  +  1 
xl=ixl 
x2=ix2 

type  Date  to  be  plotted  from  ',ixl,'  to  ix2,'  hours  ' 

call  draw_box(z,z,xpaper,ypaper) 

call  nlgpwrchg(.2) 

y=1.35 

is=8 

call  pwritxx(l. ,y, tl, is, 0,-1) 
y=y-.2 

call  pwritxx(l . ,y, t2, is, 0,-1) 
y=y-*2 

call  pvri txx(l . ,y, t8, is, 0,-1) 

y=y- • 2 

call  pvritxx(l. ,y, 1 3 , isf0,-l) 
y=y-.2 

call  pwri txx(l. ,y, t4, is, 0,-1) 
y=y--2 

call  pvritxx(l. ,y , t5, is ,0, -1) 

call  drav_box(z,z,xpaper , 1.6) 

call  pvritxx(21. , .8, t6, 30,0,0) 
xp=2 . 

ypl=xp+l . 5 

call  plotavec(xp,ypl f time, flat ,nf ' latitude' ) 


I 


I 


► 


► 


► 


» 


I 


» 


I 


» 


I 


» 


I 


Pag*  3  __LTA27 :  27-DEC-66  14:40:43  $DISK1 : ( DEHIS . PCITY . POMS 1TTPLOT . FOR ; 6 


xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl , time, f Ion, n, ' longitude' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, hdg,n, 'heading'  ) 

xp=xp+l .6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, speed, n, ' speed'  ) 

xp=xp+1.6 

ypl=xp+l .5 

call  plotavec(xp,ypl, time, flitel ,n, '441nm(log)' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, flite2,n, '488 (linear) ' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, flite3,n, 'pyrol-bulb' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, flite4,n, ' pyro2-hemi' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time,upl,n, '465nm  upwelling') 

xp=xp+1.6 

ypl=xp+l .5 

call  plotavec(xp,ypl,time,up2,n,'507nm  upwelling') 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl,time,up3,n, '532nm  upwelling') 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time,downl,n, '488nm  downwelling' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, temp, n, ' temperature' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, salin,n, 'salinity' ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl, time, fluor,n, ' fluoresence'  ) 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl , time,ptrans,n, 'X  transmission') 

xp=xp+1.6 

ypl=xp+1.5 

call  plotavec(xp,ypl , time,depth,n, 'depth' ) 


do  i=ixl,ix2 

ii=i*100 

xx=i 

write(itime, ' (i4.4) ' )  ii 

call  pvritxx(  xx,  ypl+. 15, i time, 36, 0, 0) 

enddo 

call  frame 

stop 

end 

subroutine  plotavec(yp,ypt,x,y,n, title) 
common  /paper/  xpaper,y paper, LASER 
LOGiCAL  LASER  !  true  if  ln03  —  false  if  hp7585b 
common  /xrange/  xl,x2 
character  *(*)  title 


o4 
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real*4  x(n),y(n) 
character*10  ranget , rangeb 
CHARACTER*60  STRING 

isi=36 
isil=36 

call  minmax(y,n, fmin, fmax) 
vrite(ranget, ' (f9.2)' )  fmax 
vri te(rangeb, ’ (f9.2)' )  fmin 
STRING=TITLE 
del tay=f max-f min 
deltax=x2-xl 
yuin  =deltay/(  ypt-yp  ) 
xuin  =deltax/(xpaper-2. ) 
plin  =  .2*yuin 
x2in  =  2.*xuin 
xll  =  xl+.15*xuin 
x21  =  x2-.15*xuin 

type  *, STRING 
IF  (LASER)  then 

call  set_42(l. ,xpaper-l. ,yp,ypt ,xl,x2, fmin, fmax, 1) 
isil=36 

elselF  (.NOT. LASER)  THEN 

call  set_60(l. ,xpaper-l. ,yp,ypt,xl,x2, fmin, fmax, 1) 
isil=12 

end  if 

call  draw_box(xl,  fmin,  x2,  fmax) 

call  pwritxx(xll,  fmin+plin  ,  string  ,  isil,0,-l) 

call  pvritxx(x21,  fmin+plin  ,  rangeb  ,  isil,0,  1) 

call  pwritxx(x21,  fmax-plin  ,  ranget  ,  isil,0,  1) 

call  curve(x,y,n) 
return 
end 

subroutine  curve(x,y,n) 
parameter  (hiatus  ■  20./3600.) 
real*4  x(n),y(n) 
dx=0. 
do  i=l,n 

if(i.ne.l)  dx  =  x(i)  -  x(i-l) 
xx=x(i) 

yy=y(i) 

if(i.eq.l  .or.  dx.gt. hiatus)  then 
call  frstpt(xx,yy) 

else 

call  vector(xx,yy) 

endif 

enddo 
return 
end 

subroutine  minmax(y ,n,ymin,ymax) 
real*4  y(n) 
ymin=l.e30 
ymax=-l.e30 
do  i=l,n 

if(y(i).gt.  ymax)  ymax=y(i) 
if(y(i).lt.  ymin)  ymin=y(i) 

enddo 
return 
!  end 
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c*** 
c 
c 
c 


SION 


C  AUTtH} 

C 

C  DAT® 

C 

C  REVt 
C 
C 

0****j*j***********************************************>(:*******Jlr********-**!: 

c 

parameter  (nx=80,ny=80) 
parameter  (nl=6000) 
parameter  (n2=nl*13) 
parameter  (n3=31*nl) 
parameter  (n4=14) 
parameter  (ires=25) 
parameter  (ires2=ires**2) 


PURI 


PP 


dr************************************************************** 

PROGJRAM  CONTOUR 

1)  PROGRAM  TO  CONTOUR  PANAMA  CITY  TUPS/TESS  DATA 
NORMAN  GUINASSO 
NOV.  88 
DEC.  88 


SE 

I0R 


EQUIVALENCE ( VK ( 1 ) , VK1 ( 1 ) ) 

EQUI VALENCE (VK( 1 ) , VK2( 1 , 1 )  ) 
EQUIVALENCE(ADATA( 1,1), DEPTH ( 1 ) ) 
EQUIVALENCE ( ADATA ( 1,2), TEMPER (1 ) ) 
EQUIVALENCE(ADATA( 1,3), SALIN( 1 ) ) 
EQUIVALENCE ( ADATA ( 1,4), PTRANS( 1 ) ) 
EQUIVALENCE ( ADATA ( 1 , 5) , FLUOR ( 1 ) ) 
EQUIVALENCE(ADATA( 1,6), L465NM( 1 ) ) 
EQUIVALENCE( ADATA ( 1,7), L488NM( 1 ) ) 
EQUIVALENCE ( ADATA ( 1 , 8  ) , L507NM< 1 ) ) 
EQUIVALENCE(ADATA( 1,9), L532NM( 1) ) 
EQUIVALENCE ( ADATA ( 1 , 10  ) ,TESS441( 1 ) ) 
EQUI VALENCE ( ADATA ( 1,11  ) , TESS488 ( 1 ) ) 
EQUIVALENCE ( ADATA (1,12  ) , TESSPYR 1(1)) 
EQUIVALENCE ( ADATA ( 1,13  ) , TESSPYR2 ( 1 ) ) 
EQUIVALENCE ( ADATA ( 1 , 14  ) , SECMID( 1 ) ) 


I 


1 

2 

3 


COMMON  /C0NRE1/  IOFFP,SPVAL 
COMMON  /C0NRE4/  SIZEL  ,SIZEM 
NCRT  ,ILAB 

EXT  , IOFFM 

NLM  , XLT 

common  /fopen  name/  fname 
character*60  Iname 
common  /congrid/  xmin,xmax,ymin,ymax,nnx,nny 
logical  larea(nx,ny) 


, SIZEP 
, NULBLL 
, ISOLID 
,  YBT 


,NREP, 

,  IOFFD  , 
,  NLA  , 
,  SIDE 


REAL*4  ADATA(N1,N4) 

REAL*4  FLAT(nl ) ,  FL0N(nl )  ,  DEPTH (nl^ ) 

real*4  temper(nl) , salin(nl) , ptrans(nl ) , f luor(nl) 

real*4  L465NM(nl ) , L488NM(N1 ) , L507NM(N1 ) , L532NM(N1 ) 

REAL*4  TESS441(N1 ) ,TESS488(N1 ) ,TESSPYR1(N1 ) ,TESSPYR2(N1 ) 

real*4  array(3,nl),SECMID(Nl) 

real*4  wkl(20000) , wk2(nx,ny) 

real*4  wk(n2) , iwk(n3) ,SCarr(ires2 ) 

real*4  area(nx,ny) 


logical  aver 

character*12  thedate  /'XX  June  1988'/ 
charaicter*30  title(14) 


I 


I 


> 


i 


i 


I 


> 


I 


I 


I 
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CHARACTER* 10  FMT 
character*140  in 
character*8  ctime 

conrec  common  block 
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xlt=  .11  !  xoffset 

side=. 78  !  length  of  square 

nrep=3 

nnx=nx 


nny=ny 

I0FFP=1 

SPVAL=99999. 

title(l)='DEPTH' 

TITLE ( 2 ) = ' TEMPERATURE ' 
title(3)=' SALINITY' 

TITLE(4)=' PERCENT  TRANSMISSION' 

TITLE ( 5 ) = ' FLUORESCENCE ' 

TITLE ( 6 ) = ' TUPS  465NM' 

TITLE(7)='TUPS  488NM' 

TITLE(8)='TUPS  507NM' 

TITLE(9)='TUPS  532NM' 

TITLE(10)='TESS  441NM' 

TITLE(11)='TESS  488NM' 

TITLE(12)='TESS  PYROMETER  1' 

TITLE(13)='TESS  PYROMETER  2' 

TITLE( 14)=' SECONDS' 
aver  =  .false, 
xmin  =  1000. 

ymin  =  1000 

ymax  =  -1000. 
xmax  =  -1000. 

CALL  F0PENREAD(9) 
thedate(l:2)=fname(4:5) 
type  *,  '  file  date  =',thedate 
type  '("  how  many  plots  (  max  13)  ?  '',$)' 

accept  *,nplots 
1=0 

DO  WHILE  (.TRUE.) 

M=I+1 

read(9, ' (a)' , iostat=jos)  in 
if (jos.ne.O)  goto  10 
call  pulloff(in, ctime, len) 

READ( in , * , I0STAT=I0S) 

flat(m) , f lon(m) ,depth(m) , 
temper (M) , salin(M) , ptrans(M) , f luor(M) , 
1465nm(m) , 1488nm(m) , 1507nm(m) , 1532nm(m) , 
tess441(m) , tess488(m) , tesspyrl(m) , tesspyr2(m) 


10 


if(ios.eq.O)  then 
i  =  i  +  l 

else 

type  *,  'bad  line' 

endi  f 

enddo 

continue 

k=0 

npts=i 

nret=0 

!  type  *,  npts.nl ,n4, '  CALLING  SAMEPTS  ' 

[  call  samepts( f lat , f Ion, ADATA,npts ,N1 ,N4,Nret ) 
npts=nret 

type  * , npts , '  points  left  after  samepts' 
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do  i=l,npts  I  fill  lat  Ion  arrays 
array ( 1, i)=flat(i) 
array(2, i)=flon(i) 

call  xyminmaxs (flon(i),flat(i),xmin, xmax , ymin , ymax ) 

enddo 

type  *,'  Ion  min  max  ',xmin,xmax 

type  lat  min  max  ymin, ymax 

call  mapmin(xmin, xmax, ymin, ymax) 
type  Ion  min  max  ',xmin,xmax 

type  lat  min  max  ', ymin, ymax 

ioptn=l  !  1=  mask  on  0=  mask  off 

ireach=5 

call  mask_init  (. false. , larea) 

call  mask_conrec(f lat , f Ion, npts, i reach, ioptn, larea) 

call  mask_shov( larea) 

call  pci ty_land(0, larea) 

call  mask_shov( larea) 

call’  lib$init_timer() 

DO  KN=l,nplots 

set  up  array  for  weaver: 
do  i=l,NPTS 

array(3, i)=adata(i,kn) 

enddo 

dx=xmax-xmin 
dy=ymax-ymin 
dxy=dx 

wea)v)er  parameters 

dist=  .075*dxy  !  measure  of  distance  to  fill  (  dist  must  be  <.5*djxy  ) 
ilr  =0  !  wrap  left  right 

ibt  =0  !  wrap  bottom  top 

iweav  =1  !  0=weave  lr  first;  l=veave  bt  first 
idee  =2  !  decimal  places 

ifine  =1  !  number  of  times  through  emphasize 
centr  =1.!  weight  of  center 
ichek  =0  !  i  urned  error  parameter 
wxmin  =  xmin 
wxmax  =  xmax 
vymin  =  ymin 
vymax  =  ymax 
!  nxv=nx 
nyv=ny 

type  *,  '  dist,  dx,  xy  =  '  ,dist,dx,dy 
type  *,  'nx,ny  ',nx,ny 

type  *  ,  '  min  max  ', wxmin, wxmax, vymin, wymax 

call  weaver (ar ray, np t s, ar ea, wk2 ,nxw,nyv, 
wxmin, wxmax, vymin, wymax, dist , 
ilr , ibt , iweav, idee, if ine, centr , ichek) 

type  *,  'after  weaver  ichek  =  ', ichek 

call  mask_apply(area, larea) 

call  con rec( area, nx,nx,ny ,0. ,0.,0.,0,-l,0) 


68 


i 
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j=str$triro( ti tle(kn) , title(kn) ,len) 

type  *, title(kn)(l:len) 

z0=0 . 

zl  =  xlt 

z2  =  xlt+  side 

z3  =  .05  +side 

zw=l . 

call  set  (zl,z2, .05,z3,xmin,xmax,ymin,ymax, 1) 
call  dashd( ,12,5,1) 
call  curvedh(flon,flat,npts) 
i  call  curve(f Ion, flat ,npts) 

call  pcity_land(l,larea) 
call  set  (zl,z2,zl,z2,  zl,z2,zl,z2, 1) 
call  pwri txx( . 14, . 78, title(kn) , 12, 0,-1) 
call  pwritxx( . 14, . 75, thedate, 12, 0,-1) 
iz=10 

call  map_vrite_latlon(.5,  .05-. 02  ,iz,0,0,ymin, 'lat' ) 
call  map_write_latlon( .5,  z3+.02  , iz ,0,0,ymax, ' lat' ) 

call  map_vrite_latlon(zl-.02,  .5  ,iz,90,0,xmin, 'Ion' ) 

call  map_vrite_latlon(z2+.02,  .5  ,iz,90,0,xmax, 'Ion' ) 

call  frame 

call  lib$shov_timer() 

ENDDO 

stop 

end 

subroutine  curvedh(x,y,n) 
real*4  x(n),y(n) 
logical  penup 
penup=. true. 

thresh  =  .5/60.  !  approx  .5  NM 
call  frstd(x(l),y(l)) 
do  i=2,n 

dx  =  x(i)-x(i-l) 

dy  =  y(i)-y(i-l) 

dx  =dx*cosd(y(i)) 

di  =  sqrt  (dx**2  +  dy**2) 

if  (di.gt.  thresh)  then 

type  * , ' FRSTD  n=  '  ,i 
call  lastd 

call  frstd(x(i),y(i)) 
penup  =  .false. 


endif 

enddo 

call  lastd 

return 

end 


call  vectd(x(i),y(i)) 
penup  =  .true. 


subroutine  mask_apply(area, larea) 
applys  conrec  mask  to  weaver  output  file 


!  guinasso,  circa  1988 

common  /conrel/  ioffp,spval 

common  /congrid/  xmin,xmax,ymin,ymax,nx,ny 

logical  larea(nx,ny) 

real*4  area(nx,ny) 

character*50  line 

!  i 

type  *,  '  maskapply  ',nx,ny 
do  i=l ,NX 
I  do  j=l ,NY 

if ( . not . larea( i , j ))  area( i , j )=spval 
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enddo 

enddo 

return 

end 


! 

I 


subroutine  mask  init(linit,larea) 

common  /con_gri<J/  xmin,xmax,ymin,ymax,nx,ny 

logical  larea(nx,ny) 

logical  linit 

do  j=l,ny 

do  i=l,nx 

larea(i , j )=linit 

enddo 

enddo 

return 

end 


; 


subroutine  mask  show(larea) 

common  /con_gri3/  xmin,xmax,ymin,ymax,nx,ny 

logical  larea(nx,ny) 

zl  =  .05 

z2  =  .95 

dx  =  (xmax-xmin)/float(nx-l) 

dy  =  (ymax-ymin)/float(ny-l) 

call  set(zl, z2,zl, z2,xmin,xmax,ymin,ymax, 1) 

do  j=l,ny 

ynode=ymi n+ ( j - 1 ) *dy 
do  i=l,nx 

xnode=xmin+(i-l)*dx 

if(larea(i, j))  call  ticm(xnode,ynode,dx/15. ) 

enddo 

enddo 

call  frame 
end 


subroutine  mask_conrec(flat , f Ion, npts, i reach, ioptn,larea) 
common  /congrid/  xmin,xmax,ymin,ymax,nx,ny 
logical  larea(nx.ny) 
flat  flon  defines  a  track 
ioptn=  1  turns  on  grid  pts  near  track 
ioptn=  0  turns  off  grid  points  near  track 
xmin,xmax,nx  ymin,ymax,ny  defines  a  grid  nx  by  ny 

lflag(i,j)  is  true  if  point  i,j  is  within  distl  of  any  point  on  track 
this  defines  a  a  mask  that  can  be  applied  to  a  data  grid  passed  to  conlrec 
guinasso,  circa  1988 


real*4  flat(npts) , flon(npts) 
logical  maskl,mask2 
dist(x,y)=  (x**2  +  y**2) 

if ( ioptn.eq. 1)  then 
maskl=  .true. 
inc=  1 

else 

maskl=  .false, 
inc  =  -1 

endi  f 

)  itot=nx*ny 
tot=i tot 

dx=(xmax-xmin)/f loat(nx-l) 
dy=(ymax-ymin)/f loat(ny-l) 
j  |  dis  =  min(dx.dy) 
distl  =  (ireach*dis)**2 


!  rem 
!  N  « 
!  MD 
!  ND 
!  gui 
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TYPE  *, 'ENTERING  MASKCONREC  with  ioptn  =  ',ioptn,inc 
type  '<"  dx,dy  distl="  ,3fl2.4)' ,  dx,dy,distl 
type  *, '  ' 
ion=0 
do  j=l,ny 

ynode=ymin+( j-l)*dy 
do  i=l,nx 

xnode=xmin+(i-l)*dx 

i f (maskl . eq . larea( i , j ) )  then 

else 

do  k=l,npts 

dely  =  flat(k)-ynode 
delx  =  flon(k)-xnode 
dist2  =  dist(delx,dely) 
if(dist2  .It.  distl)  then 
larea(i,j)  =  maskl 
ion  =  ion  +  inc 
goto  5 

else 
end  if 

enddo 
end  if 

continue 

f ion=100.*ion/tot 

if(mod(i,ny).eq.O  .and.  mod(j ,10) .eq.O)  type  109,j,fion 

enddo 

enddo 

type  *,  '  Mask  conrec-i-  cells  changed,  total  cells=  ',ion,itot 
return 

format('+  j,  percent  switched  =  ',15,f6.1,'  ') 

end 

subroutine  samepts(flat, flon,ADATA,n,ndl,nd2,nret) 

)ves  data  with  duplicate  locations 
NUMBER  OF  LINES 
=  DIMENSION  OF  ADATA 
■  NUMBER  OF  DATA  VARIABLES 
lasso,  circa  1988 

real*4  f lat(n) , flon(n) ,ADATA(ndl .nd2) 
logical  skip(6000) 
ireject=0 
do  i=2,n 

do  j=l , i-1 

if (flat(i).eq. flat(j)  .and.  flon(i) .eq. flon( j ))  then 
ireject  =  ireject+1 
skip(i)  =  .true, 
goto  10 


end  if 


enddo 


enddo 

k=0 

do  i=l,n 


skip(i )=. false. 


IF(SKIP(I) )  THEN 
else 

f lat(k+l)=f lat( i ) 
f lon(k+l)=f lon( i ) 

DO  KD=1 ,ND2 

ADATA (k+1 ,KD) = ADATA ( i , KD ) 


endi  f 


ENDDO 

k=k+l 

nret=k 
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enddo 

type  total,  returned,  rejected  ,n,nret,ireject 

return 

end 

subroutine  pci ty_land(iflag, larea) 
common  /con_grid/  xmin,xmax,ymin,ymax,nx,ny 
logical  larea(nx,ny) 
real*4  f la t ( 500 ) , f Ion ( 500 ) 

open(19,name='pcitymap. flO' , status®' old' ,iostat=kos) 

if(kos.ne.O)  stop  ' pcitymap. f 10  not  found' 

ipold=0  !  last  line  # 

npts=0 

nlines=0  !  lines  plotted 

do  while  (.true.) 

read ( 19 , * , ios  ta t=kos )  num , i p , x , y , alAt , alon 
alon=-alon 

if (ip.ne. ipold  .or.  kos.eq.-l)  then 
if (nlines.ne.0)  then 

if (iflag.eq.l)  then 

call  tcurve(flon,flat,npts)  - 
else 

call  mask_conrec(flat,f Ion, npts, 1,0, larea) 

subroutine  mask_conrec( flat , f Ion, npts, ireach, ioptn, larea) 

end  if 

if (kos.eq.-l)  then 
close  (19) 
return 

endif 

endif 

nlines=nlines+l 

npts=l 

flat(l)=alat 

flon(l)=alon 

ipold=ip 

type  *,  '  pcity_land-i-nlines  =  ',nlines 

else 

npts=npts+l 
f lat(npts)=alat 
f lon(npts)=alon 

endif 

enddo 

end 

subroutine  mapmin( f lonl , f lon2, flatl,flat2) 

if(flonl.gt.O)  stop  'mapmin' 

flonl=-flonl 

f lon2=-f lon2 

iflatl  =  flatl*60. 

iflat2  =  flat2*60. 

iflonl  =  flonl*60. 

iflon2  =  flon2*60. 

f latl=f loat( if latl  )/60. 

f lat2= float (if lat2+l)/60. 

f lonl =- float ( if lonl +l)/60. 

flon2=-float(iflon2  )/60. 

return 

end 

subroutine  tcurve(x,y ,n) 
external  between 

common  /congrid/  xmin,xmax,ymin,ymax,nx,ny 
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real*4  x(n),y(n) 

LOGICAL  between 
logical  penup 
penup  =  .true, 
do  i=l,n 


5 


if(  betveen(xmin,xmax,  x(i))  .and. 
betveen(ymin,ymax,  y(i)) 

)  then 

if(penup)  then 

call  frstpt  (  x(i),y(i) 
penup  =  .false. 

else 

call  vector  (  x(i),y(i) 

end  if 

else 

penup  =.true. 

end  if 


) 

) 


enddo 

return 

end 
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2****^*************************************************************** 


1)  sets  up  map  for  land  boundary  plotting 


NORMAN  GUINASSO 


circa  1985 


REVISION 


NOV.  1988 


2S  calibration  data  and  x,y  station  data  from  a  *.mca  file 
Ls  mapcalib  to  make  a  *.cof  file  from  calibration  data 
)  can  call  map_make_square  to  make  a  *.cof  file  for  a  square  projection 
3  looks  for  a  *.flt  auxiliary  file  for  data  input  as  lat  Ion 
i  calculates  x  y  or  lat  long  for  each  data  point 
vrites  all  to  a  *.flO  file 

files  take  their  *  name  from  the  specified  input  file 

k ******************************************************************* 


common  /mapname/  infile 
common  /mapof fsets/  xoff ,yof f , fac 
character*70  title 
character  ans 
!  calibration  data: 

real*4  x(110) ,y(110)  ! 

real*4  flat(HO) , flon(llO)  ! 

!  data  points: 

integer*4  ident(900) 
real*4  xx(900) ,yy(900) ,xz,yz  ! 

real*4  xl(900) ,yl(900)  ! 

logical  square  ! 

logical  finnmark  /.FALSE./  ! 

logical  digit_vith_calib 
logical  toxy  /.true./ 

logical  to  latlon  /.false./ 


!  x,y  coordinates  in  inches 
!  lat  and  long  in  decimal  degree 


!  x,y  for  data  points 
!  long,  lat  fro  data  points 
!  if  square  plot  set  up  a  simple  map 
!  special  handling  for  finnmark 


character*30  inf ile, out file, auxfile 
character*60  auxformat 
logical  rotelkwise 
logical  fit  file 


calibrated  or  square 
x,y  points  or  fit  file 

|  fac=l. 

j  j  =  1 i b$ ini t_timer( ) 

type  *,  'Program  Ml  version  1.1' 


square=. true. 

type  '("  calibrated  or  square?  (c/s)  — >",$)' 
accept  '(a)',  ans 

if(ans.eq. 'c'  .or.  ans.eq.'C')  square=. false. 

I  if(square)  then 

type  *,'  scale?  ' 

!  accept  *,sqscale 

!  end if 


format ( '  Enter  input  file  name  ->  ',$) 
type  101 

accept  103, infile 


75 
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II 


call  chgext(inf ile,outf ile, ' £10' ) 
open(10,name=outf ile, err=991, status* 'NEW' , 
carriagecontrol='LIST' ) 


! - 

!  if 
!  f  iri 
f  wi 
!  stal 


tW 


fit  file  is  found,  then  assume  lat  Ion  is  there 
t  record  contains  format 
each  of  the  following  records  containing: 
ion  lat  deg  lat  min  flat  secs  Ion  deg  Ion  min  flon  secs 


call  chgext(infile,auxfile, '  fit' ) 
open(ll,name=auxf ile, status* 'OLD' ,err=10) 

!  foujnjd  aux  file: 

digit_vith_calib  =  .false. 

type  *,  'Found  auxiliary  file  with  lat  Ion  data:  '  ,auxfi 
read(ll, ' (a)'  )  auxformat 

type  '(lx, a)','  ',  auxformat 

goto  11 
|not  find  aux  file: 

digit_with^calib  *  .true. 


le 


!  didj 
10 
! 

11 
!  — 

!  stal: 


continue 


t: 


fltfile*  .not.  digitwithcalib 

if ( .not. (flt_file. and. square))  then 

we  need  to  read  calibration  data 

call  chgext(infile, infile, 'mca' ) 

OPEN (9 , name=inf i le , er r*990 , s  ta tus* ' old ' , readonly ) 
read(9,100)  n, title 

type  109,  n, title 

write(10, 109)  n, title 

if  (n.gt.110)  then 

type  *,  '  this  will  no l  work,  n  too  big' 
stop 

endif 


!  rea|d| 
138 


! 


in  calibration  data 

format( '  Rotate  clockwise?  (y/n)  — >  ',$) 
type  138 

accept  '(a)',  ans 
rotelkwise  =  .false. 

if (ans.eq. 'y'  .or.  ans.eq.'Y')  rot_clkwise=. true, 
do  i=l,n 

read(9,*,end=900,err=901)  Lat,Long,x(i),y(i) 
if (rotelkwise)  then 
xtemp  =  x(i) 
x(i)  =  y(i ) 
y(i)  =  -xtemp 

endif 

TYPE  *, 'DENIS  MODIFIED  ME' 
flat(i)  =  ghexij(lat) 
flon(i)  =  ghexij(long) 

TYPE  102, flat(i),flon(i),x(i),y(i) 
wri te(10, 102)  f lat(i ) , f lon(i) ,x(i) ,y(i ) 

enddo 

endif 


I 

I 


a-a. 
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if  (square)  then 

call  ioap_niake_square(sqscale,  1)  I  (scale,  1=  write  file)) 

else 

call  map_calib(flat,flon,x,y,n) 

end  if 


if(finnmark)  then 
nskip=21 
do  k=l,nskip 
read(9,103)  title 
enddo 

end  if 


i=l 

npoints=0 


!  next  point  to  read 
!  points  read 


if (digit_with_calib)  then 

read  (9  ,103)  title 
wri te(10, 106)  title 

read  station  locations  from  digitizer  file  as  x,y 
do  while  ( . true. ) 

read(9,*,end=50,err=50)  ident(i) ,xx(i) ,yy(i) 
if (rotclkwise)  then 

xtemp  =  xx(i) 

xx(i)  =  yy(i) 

yy(i)  =  -xtemp 

end  if 

type  107 , ident(i) ,xx(i) ,yy(i) 
if (finnmark)  ident(i)=ident(i)+2000 
npoints  =npoints+l 
xt=xx(i) 
yt=yy(i) 

call  map_xform(ft,fg,xt,yt,to_latlon)  !  to  lat 
xl(i)=fg 
yi(i )=f t 

write(10, 108)  i , ident(i) ,xx(i) ,yy(i) ,yl(i ) ,xl(i ) 
i  =  i+1 

enddo 

else 

read  station  locations  from  auxiliary  file  as  lat  long 
do  while  ( . true. ) 

read ( 1 1 , auxf orma t , end=50 , err  =44 )  iden  t ( i ) , 
latdeg, latmin,secslat , londeg, lonmin, secs Ion 
npoints  =npoints+l 

type  '(2i4,f6.1,2i4,f6.1)', 
latdeg, latmin,secslat , londeg, lonmin, secslon 


long 


76 


funjcjtion  ghexijf  converts  degrees  minutes  and  seconds  to  floating  degrees 
degrees  and  minutes  are  integers,  seconds  are  floating  point  ' 

fg=ghexijf (londeg, lonmin, secslon) 
f t=ghexi j  f ( latdeg, lat min, secs lat) 
type  *,ft,fg 

type  ' ( i4, 2f 10. 3) ' ,  iden t ( i ) , f t , fg, i den t ( i ) , xz ,yz  1 

call  map_xform( f t , fg,xz ,yz , to_xy )  !  to  x,y  I 

xl(i)=fg  !  long 

yl(i)=ft  !  lat  1 

xx( i  )=xz 
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» 


» 


yy(i)=yz 

write(10,108)  i,ident(i),xx(i),yy(i),yl(i),xl(i) 
i  =  i+1 


44 

45 


goto  45 

type  *,  '  error  while  reading  on  unit  11  ',auxfile 
enddo 

endif!  (digit  withcalib) 


» 


» 


900 

901 

990 

991 

100 

102 

103 

106 

108 

107 

109 

!  squa 
!  squa 
!  cali 
!  cali 


continue 

type  *,  'number  of  points  =  '.npoints 

j=lib$shov_timer() 

stop  'normal  ending' 

stop  'program  ml-f-unexpected  end  of  data' 
stop  'program  ml-f-error  reading  calibration  data' 
stop  'program  ml-f-error  opening  input  file  unit  9' 
stop  'program  ml-f-error  opening  output  file' 

format  (i3,a)  !  read  format 

format(lx,2fl0.3,2fl0.4) 

format(a) 

format(lx,a) 

format(lx,2I7,4fl2.5) 

format (lx, i4,2f8.2) 

format(lx, i 3 , '  - >',a) 

re  with  calibration  data  and  point  data 
re  without  calibration  data  and  with  fit  file 
bration  data  and  x,y  data 
bration  data  and  lat  long  data  in  fit  file 
end 


I 


78 


I 


I 


I 


> 


> 


) 
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/*  PRi 
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m 


PUR] 


/* 

/* 

/* 

/* 

/* 

/*  DA) 


RAM 
POSE 
AUtTHOR 


A*********************************************************/ 

*/ 

READ . SAS 


We 


1)  sas  program  to  consolidate  .fin  files 
NORMAN  GUINASSO 
NOV.  88 
DEC.  88 


I/* 

/*  REVISION 
/* 

/****)*j********************************  ★★★★**********★***★**★★**★*  / 


*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 


FILENAME  PCITY1  ' [DENIS. PCITY] 19 JUN88 . FIN' ; 

LIBNAME  NORDA  ' [DENIS. PCITY] ' ; 

DATA  tempi; 

INFILE  PCITY1; 

INPUT  YEAR  MONTH  DAY  TIME  $  SECSMID 
LAT  LON  T_DELAY1  T_DELAY2 
HEADING  SPEED  VREF  TESS441  TES448 
TESSPYR1  TESSPYR2  TEMPER  SALIN 
PTRANS  FLUOR  L465NM  L507NM  L532NM 
TRANSVLT  FL_SIG_V  FL_SCL_V  L488NM  DEPTH; 

data  temp2; 

FILENAME  PCITY2  ' [DENIS. PCITY] 20 JUN88. FIN' ; 

INFILE  PCITY2; 

INPUT  YEAR  MONTH  DAY  TIME  $  SECS_MID 
LAT  LON  TDELAY1  T_DELAY2 
HEADING  SPEED  VREF  TESS441  TES448 
TESSPYR1  TESSPYR2  TEMPER  SALIN 
PTRANS  FLUOR  L465NM  L507NM  L532NM 
TRANSVLT  FL_SIG_V  FL_SCL_V  L488NM  DEPTH; 

data  temp3; 

FILENAME  PCITY3  ' [DENIS. PCITY] 21 JUN88. FIN' ; 

INFILE  PCITY3; 

INPUT  YEAR  MONTH  DAY  TIME  $  SECS_MID 
LAT  LON  T  DELAY1  T_DELAY2 
HEADING  SPEED  VREF  TESS441  TES448 
TESSPYR1  TESSPYR2  TEMPER  SALIN 
PTRANS  FLUOR  L465NM  L507NM  L532NM 
TRANSVLT  FL_SIG_V  FL_SCL_V  L488NM  DEPTH; 

data  norda.all; 

j  set  tempi  temp2  temp3; 

I  PROC  MEANS; 

ENDSAS; 


i 


l 


l 


l 
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/*** 
/* 

/*  PRl 


PUW 


POSE 


/* 

V* 

/* 

/* 

/* 

/* 

/*  AUTHOR 
/* 

/*  DAtl 


/* 

/* 

/* 

/*** 


OiGRAM 


WE 


REM: 


**********************************************************/ 

*/ 

PUTALL.SAS 


I  SION 


1)  sas  program  to  make  smooth  data  file 

2)  scans  large  sas  file  and  outputs 
selective  records 

NORMAN  GUINASSO 

NOV.  88 

DEC.  88 


*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 


r ***************************** ********★********************/ 

FILENAME  PCITY  ' [DENIS. PCITY]all21.nav ' ; 

LIBNAME  NORDA  ' [DENIS. PCITY] ' ; 

DATA  TEMP; 

FILE  PCITY; 

SET  NORDA. all; 

IF  DAY=  21  THEN  PUT 

TIME  1-10  LAT  10.6  LON  11.6  DEPTH  6.2 

@40  TEMPER  SALIN  PTRANS  FLUOR  L465NM  L488NM  L507NM  U532NM 
TESS441  TES448  TESSPYR1  TESSPYR2  SECSMID; 

ENDSAS; 
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I 
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„  __ ^ .  •  • _ • _ _ __  _  _  •  _ ... •  ____ •  _ 

•  —  a  •  • - • - •  —  —  —  #  •  •  •  •  • 

It**********************************************************/ 

*/ 

3RAM  CORR. SAS  */ 

*/ 

POSE  1)  sas  program  to  produce  correlation  matrix  */ 

*/ 

10R  NORMAN  GUINASSO  */ 

*/ 

3  DEC.  88  */ 


fllSION 


DEC.  88 


t ***************************** ***************************** / 

LIBNAME  NORDA  ' [DENIS. PCITY] ' ; 
data  temp; 

set  norda.all; 
tess448  =  tes448; 

keep  tess448  tess441  tesspyrl  tesspyr2 
temper  salin  ptrans  fluor 
L465NM  L507NM  L532NM  L488NM  DEPTH; 
proc  corr  data=TEMP; 

ENDSAS; 


fl'ftffrtonrtonrtoortft 
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|C*** 
progra 

PURPO 


*w **************************************************************** 


c 


*W******************************************************************f 


c 

c  opt] 
iC 


M  READ  350 

SE  1)  RT-11  F0RTRAN_IV  program  to  read  disk  to  recover  filejs 

2)  modified  for  NORDA  PRO-350  1988 
AUTHOR  NORMAN  GUINASSO 

DATE  1982 

REVISION  OCT.  1988 


integer*2  dblk(4) 
integer*2  dblko(4) 
byte  buff(514) 
byte  bell 
byte  ans 
logical  octal 
logical  ascii 
logical  vtf 

data  bell  111 
data  dblk  /3RDU0.3R 


i 

! 

,3R 


floppy 
output  file 
block  buffer 
ring 

a  character 

print  to  screen  in  octal 
print  to  screen  in  ascii 
write  to  file 


,3R 


data  dblko/3RDZl , 3RREA, 3RDXX , 3R0UT/ 


!  the  pro350 
!  duO:readxx.out 


buf f (513)=0 
buf f (514)=0 

type  *,  '  program  to  use  hardware  read  to  recover  files' 
type  *, 


lb 


bns: 


ascii 

wtf 

type  *, 


from  floppy  -  -  -  guinasso,  circa  1987' 

=  1+istop-istart  !  length  of  output  file 

=  .false.  !  write  first  75  characters  of  block  asclii 

=  .false,  f  write  first  75  characters  of  block  asqii 

'  printout  block  in  Ascii,  Octal,  or  None?' 
accept  101, ans 


'A')  ascii 
'0')  octal 


if(ans.eq. 'a'  .or.  ans.eq 
if (ans. eq. 'o'  .or.  ans.eq 
type  *,'  write  to  file?  ' 
accept  101, ans 

if (ans.eq. 'y'  .or.  ans.eq. 'Y')  wtf 


=. true . 
=. true. 


= . true. 


!  channel  for  floppy 
!  channel  for  output  file 


ic  get  channels 

ich=igetc( ) 

i  ichout=igetc( ) 

[c  1 

|c  opein  disk  as  file: 

I  !|  if ( ifetch(dblk) . It .0)  stop  'bad  fetch'  !  fetch  floppy  handler 
[  if(lookup(ich,dblk).lt.0)  stop  'lookup  failed'  !  open  floppy 

c 

|c  if(.not.wtf)  goto  4  !  writing  to  file 

i  i f ( ienter(ichout ,dblko, lb) . It .0) 

%  stop  'enter  failed'  !  open  output 

type  file  opened  for  output' 

iblko=0  !  first  block  is  0 

c 


O  CO 
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istart=0 

istop=11000 

iblk=istart 


14*03:10  $DISK1 : ( DENIS . PCITY ] READS 50 . FOR; 3 

_ * _ _  _ _ * _ _ _ •  _ . _ * 


£ 


|C 

|c  harjdl 


do  10  i=l, 10100 
ierr=0 
continue 

.  •  j=ireadw(256, buf f , iblk, ich) 

if(j.eq.-l)  goto  11  !  end  of  file 

if(j.eq.-2)  goto  8  !  error 

type  •  BLOCK  =  iblk 

iblk  =  iblk+1 

jin=index(buff , '20-JUN-88' ) 
wtf=jin 

if(wtf)  j  =  ivritv(256,buff,iblko,ichout) 

if(wtf)  iblko  =  iblko+1 

if(octal)  type  100,(buff(k),k=l,75) 

if(ascii)  type  103,(buff(k),k=l,75) 

if(iblk.gt.istop)  goto  20 
goto  10 

error  handling 


10 

11 

20 


100 

101 

102 

103 


type  102  , bell, iblk 
ierr=ierr+l 

if (ierr.lt. 4)  goto  5  !  retry  3  times 

stop  '  hardware  error' 


continue 
goto  20 

type  *,'eof,  block=' ,  iblk 
call  iclose(ich) 
call  iclose(ichout) 
stop 

forma t(20o4) 
format(a) 

format(lx,a, '  hard  error,  block=',i4) 

format(80al) 

end 
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It**********************************************************/ 

*/ 

SRAM  PLOT.  SAS  */ 

*/ 

POSE  1)  sas  program  to  plot  light  data  */ 

*/ 

10R  NORMAN  GUINASSO  */ 

*/ 

Z  DEC.  88  */ 

*/ 

/ilSION  DEC.  88  */ 

*/ 

********************* ********** **************************** / 
r  LIGHT  DATA  AND  DEPTH*/ 

LIBNAME  NORDA  ' [DENIS. PCITY] ' ; 

DATA  NORDA. ALLSUB; 

set  norda.all; 
tess448  =  tes448; 

NORMUP  =  L507NM/TESSPYR1 ; /*NORMALIZED  UPWELLING  L507*/ 
INORMUP  =  l./N0RM_UP; 

KEEP 

,  DAY  tess448  tess441  tesspyrl  tesspyr2 

temper  salin  ptrans  fluor  N0RM_UP  INORMJUP  - 
L465NM  L507NM  L532NM  L488NM  DEPTH; 

DATA  TEMP; 

SET  NORDA. ALLSUB; 

IF  DAY=19  AND  DEPTH<12.; 

PROC  PLOT  DATA=TEMP ; 

PLOT  DEPTH*INORM_UP ; 

TITLE  'INVERSE  NORMALIZED  507NM  UPWELLING  IRRADIANCE' ; 
TITLE2  '19  JUNE  1988'; 

DATA  TEMP; 

SET  NORDA. ALLSUB; 

IF  DAY=20  AND  DEPTH<12.; 

PROC  PLOT  DATA=TEMP ; 

PLOT  DEPTH*INORM_UP; 

TITLE  'INVERSE  NORMALIZED  507NM  UPWELLING  IRRADIANCE'; 
TITLE2  '20  JUNE  1988'; 


DATA  TEMP; 

SET  NORDA. ALLSUB; 

IF  DAY=21  AND  DEPTH<12.; 

PROC  PLOT  DATA^TEMP; 

PLOT  DEPTH*INORM_UP ; 

TITLE  'INVERSE  NORMALIZED  507NM  UPWELLING  IRRADIANCE'; 
TITLE2  '21  JUNE  1988'; 


ENDSAS; 


APPENDIX  B 


CONTOUR  MAPS 


DEPTH 


85 


M.Ofr.SS 


M.9t,S8 


3(r02'N 

CONTOUR  FROM  2.0000  TO  21.000  CONTOUR  INTERVRL  OF  1.0000  PT (3. 3) =  99999. 


30“06'N 


86 


*,o*.se 


A.9*  ,90 


30R02'N 

CONTOUR  FROM  26.900  TO  27.800  CONTOUR  INTERVAL  OF  0.50000E-0!  PT (3. 3) =  99999. 


SALINITY 


87 


A.0t.S9 


H,9V  0S8 


30*02'N 

CONTOUR  FROM  33.870  TO  3-1.100  CONTOUR  INTERVRL  OF  0. 10000E-01  PT (3. 31  =  99999. 


PERCENT  TRANSMISSION 


88 


M,0fr<,S8 


*.9nS8 


3Cf02’N 

CONTOUR  FROM  86.000  TO  92.400  CONTOUR  INTERVPL  OF  0.40000  PT ( 3. 3) =  99999. 


FLUORESCENCE 


89 


fc.0t-oS8 


n.gt’oSg 


3(T02'N 

CONTOUR  FROM  98.000  TO  210.00  CONTOUR  INTERVAL  OF  7.0000  PT(3.3I=  99999. 


465NM 


30"02'N 

CONTOUR  FROM  I. OOOO  TO  25.000  CONTOUR  INTERVAL  OF  1.0000  PT13.31 


488NM 


91 


30°  02'  N 

CONTOUR  FROM  12.900  TO  18.900  CONTOUR  INTERVAL  OF  0.30000  PT(3. 31 =  99999. 


TUPS  507NM 


92 


3(f  02'N 

CONTOUR  from  1.0000  TO  33.000  CONTOUR  1NTERVRI  OF  1.0000  PT C 3. 3) =  99999. 


TESS  441NM 


k.SPSQ 


30"02'N 

CONTOUR  FROM  12.800  TO  16.800  CONTOUR  INTERVAL  OF  0.20000  PT13.3)-  99999. 


TESS  488NM 


30“02'N 

CONTOUR  FROM  16.670  TO  16.708  CONTOUR  INTERVAL  OF  0.20000E-02  PT 13. 31 


TESS  PYROMETER  2 


i 


DEPTH 


I 


30°03'N 

CONTOUR  FROM  3.0000  TO  21.000  CONTOUR  INTERVAL  OF  1.0000  PT(3.3) 


30*08'N 


30*03'N 

CONTOUR  FROM  85.600  TO  93.200  CONTOUR  INTERVAL  OF  0.10000  PT ( 3. 3) =  99999. 


30*08'N 


I 


101 


30-03'N 

CONTOUR  FROM  33.890  TO  3*3.130  CONTOUR  INTERVAL  OF  0.  10000E-01  PT (3.31=  99999. 


102 


I 


3(T03'N 

CONTOUR  FROM  102.00  TO  153.00  CONTOUR  INTERVAL  OF  3.0000  PT ( 3. 31 =  99999. 


30“  08' N 


I 


103 


M.St  .90 


30"03'N 

CONTOUR  FROM  1.0000  TO  21.000  CONTOUR  INTERVAL  OF  1.0000  PT ( 3. 31 =  99999. 


30“  08'N 


i 


30CO3N 

CONTOUR  FROM  15.500  *'  18.400  CONTOUR  I NTERVRL  OF  0.10000  PT (3. 3) =  99999. 


30*  08’N 


105 


30*03'N 

CONTOUR  FROM  2.0000  TO  27.000  CONTOUR  INTERVAL  OF  1.0000  PT (3. 31 = 


30*08'N 


106 


I 


3(T03'N 

CONTOUR  FROM  1.0000  TO  28.000  CONTOUR  INTERVAL  OF  1.0000  PT (3. 31 =  99939. 


I 


30°  08'N 


3(f03'N 

CONTOUR  FROM  16.677  TO  16.728  CONTOUR  INTERVAL  OF  0.30000E-02  PT ( 3. 3) =  99999. 


30‘08'N 


§ 


30'03'N 

CONTOUR  FROM  7.3800  TO  8.8200  CONTOUR  INTERVAL  OF  0.90000E-01  PT(3.3)=  99999. 


METER  2 


I 


3<F03'N 

CONTOUR  FROM  7.3600  TO  8.8000  CONTOUR  INTERVAL  OF  0.80000E 


DEPTH 


111 


30“05’N 

CONTOUR  FROM  3.5000  TO  16.100  CONTOUR  INTERVRL  OF  0.70000  PT(3.3!=  99999. 


30"  09'N 


112 


I 


SO’OS'N 

CONTOUR  FROM  27.600  TO  29.010  CONTOUR  INTERVAL  OF  0.80000E-01  PTI3.3)-  99999. 


30“  09'N 


30“05'N 

CONTOUR  FROM  31.500  TO  33.900  CONTOUR  INTERVAL  OF  0.10000  PT(3. 31 =  99999. 


30*  09'N 


3(r05N 

CONTOUR  FROM  78.100  TO  90.300  CONTOUR  INTERVAL  OF  0.70000  PT ( 3. 31 =  99999. 


30s  09'  N 


115 


I 


30*05’N 

CONTOUR  FROM  It 0.00  TO  S40.00  CONTOUR  INTERVRL  OF  20.000  PT ( 3. 31 =  99999. 


30*09'N 


1  16 


*.6E  .96 


M.St’.Se 


3(f  05N 

CONTOUR  from  O.OOOOOE>00  TO  6.8000  CONTOUR  INTERVAL  OF  0.40000  PT ( 3. 31 =  99999. 


30*09'N 


117 


3<r05'N 

CONTOUR  FROM  7.5000  TO  15.500  CONTOUR  INTERVAL  OF  0.50000  PT(3.3t=  99999. 


3Cf  09'N 


30"05'N 

CONTOUR  FROM  O.OOOOOE+OO  TO  10.200  CONTOUR  INTERVRL  OF  0.60000  PT(3.3)=  99999. 


30“  09'N 


119 


M.6C  .88 


\ 


30*  05'N 

CONTOUR  FROM  O.OOOOOE*QO  TO  10.200  CONTOUR  INTERVRL  OF  0.60000  PT(3.3I-  99999. 


3<f09'N 


120 


3<r05'N 

CONTOUR  FROM  8.4000  TO  16.000  CONTOUR  INTERVAL  OF  0.40000  PT(3.3)=  99999 


30*09'N 


30*05'N 

CONTOUR  FROM  16.638  TO  16.698  CONTOUR  INTERVPL  OF  0.30000E-02  PT ( 3. 3) =  99999. 


3<y*09'N 


3CT  05N 

CONTOUR  FROM  1.2000  TO  8.4000  CONTOUR  INTERVAL  OF  0.40000  PT(3.3I=  99999. 


30“  09' N 


30"05'N 

CONTOUR  FROM  1.2000  TO  8.8000  CONTOUR  INTERVAL  OF  0.40000  PT(3,3)=  99999. 


t 

124 


> 


> 


> 

APPENDIX  C 

» 


VAX  9-TRACK  TAPE  LISTING 


» 


i 


arnone.lis 


Thu,  Dec  28, 1989 


1 


Listing  of  save  set(s) 


Save  set : 

Written  by: 

UIC : 

Date : 

Command: 

Operating  system: 
BACKUP  version: 
CPU  ID  register: 
Node  name: 

Written  on: 

Block  size: 

Group  size: 

Buffer  count: 


PCITY. BCK 
DF.NIS 

[000200,000007] 

28-DEC-1988  14:30:50.31 

BACKUP  [DENIS. PCITY. ..] /LIST=ARNONE. LIS  MUB0 : PCITY . BCK/SAVE_SET 
VAX/VMS  version  V4 . 7 
V4.7 
08000000 
_GERGA : : 

_GERGA$MUB0 : 

8192 

10 

3 


[DENIS. PCITY] 19J1600.CNV; 3 
[DENIS. PCITY] 19J1600.TIM;2 
[DENIS . PCITY] 19J1652 . CNV; 3 
[DENIS. PCITY] 19J2022.CNV; 3 
[  DENIS . PCITY ] 1 9JUN88 . FIN; 1 
[DENIS. PCITY] 19JUN88.PTO;l 
[ DENI S . PCITY ] 1 9JUN8 8 . TES ; 1 
[  DENI S . PCI TY ] 1 9JUN8 8 . TUP ; 3 
[DENIS. PCITY] 19JUN881 . RDJ;  1 
[DENIS. PCITY] 19 JUN882 . RDJ; 1 
[DENIS . PCITY] 20 J1505 . CNV;  1 
[DENIS . PCITY] 20 JUN88 . FIN; 2 
[DENIS. PCITY] 20JUN88 .PTO; 1 
(DENIS . PCITY] 20 JUN88 . TES;  1 
[DENIS. PCITY] 20JUN88. TUP; 1 
[DENIS . PCITY] 20 JUN881 . RDJ;  1 
(DENIS . PCITY] 21 J1155 . CNV;  1 
[DENIS. PCITY] 21J1250.CNV;1 
[DENIS. PCITY] 21 J1353 . CNV; 1 
[DENIS. PCITY] 21 J1441. CNV; 1 
[DENIS. PCITY] 21 J1519 . CNV;  1 
[  DENIS . PCITY ] 21 J1 602 . CNV; 1 
[DENIS . PCITY] 21 J1630 . CNV; 1 
[DENIS. PCITY] 21 J1742 . CNV; 1 
[DENIS. PCITY] 21J1839.CNV;1 
[DENIS. PCITY] 21 JUN88 . FIN;  1 
[DENIS . PCITY] 21 JUN88 . PTO; 1 
[DENIS . PCITY] 21 JUN88 . TES;  3 
[DENIS. PCITY] 21 JUN88 . TUP; 1 
[DENIS. PCITY] 21JUN881 .RDJ; 1 
[DENIS . PCITY] 21 JUN882 . RDJ; 1 
[DENIS. PCITY] 21 JUN883 . RDJ;  1 
[  DENIS. PCITY] 21 JUN88 4 .RDJ; 1 
[  DENIS. PCITY] 21 JUN885. RDJ; 1 
[DENIS.PCITYJADDNUM.EXE;  3 
[DENIS. PCITY] ADDNUM. LIS; 4 
[  DENIS . PCITY ] ALL . SSD ; 2 
[  DENIS . PCITY ] ALL1 9 . NAV; 1 
[  DENI S . PCITY ] ALL20 . NAV ; 1 
[DENIS . PCITY ] ALL21 . NAV; 1 
[DENIS. PCITY] ALLSUB. SSD;  2 
[DENIS. PCITY] ARNONE.LIS;  1 
[DENIS . PCITY ] BEND . OUT;  59 
[DENIS . PCITY ] BEND . OUT;  58 


166 

12- JUL-1988 

12:09 

31 

29-JUN-1988 

08:41 

655 

12-JUL-1988 

12:13 

287- 

12-JUL-1988 

12:14 

1727 

28-NOV-1988 

11:41 

933 

28-NOV-1988 

11:39 

565 

ll-JUL-1988 

11:00 

1106 

12-JUL-1988 

12:34 

386 

ll-JUL-1988 

09:51 

179 

ll-JUL-1988 

09:54 

226 

3-NOV-1988 

09:19 

353 

28-NOV-1988 

11:18 

191 

28-NOV-1988 

11:17 

247 

ll-JUL-1988 

11:05 

226 

3-NOV-1988 

09:26 

247 

ll-JUL-1988 

10:01 

158 

12-JUL-1988 

12:39 

196 

12- JUL-1988 

12:40 

118 

12-JUL-1988 

12:40 

63 

12-JUL-1988 

12:41 

106 

12-JUL-1988 

12:41 

61 

12-JUL-1988 

12:42 

170 

12-JUL-1988 

12:42 

63 

12-JUL-1988 

12:43 

149 

12-JUL-1988 

12:43 

1681 

28-NOV-1988 

12:01 

908 

28-NOV-1988 

11:59 

577 

18-JUL-1988 

14:03 

1077 

12-JUL-1988 

12:55 

112 

ll-JUL-1988 

10:02 

77 

ll-JUL-1988 

10:03 

136 

ll-JUL-1988 

10:04 

140 

18-JUL-1988 

13:33 

116 

18- JUL-1988 

13:35 

6 

2-DEC-l 988 

08:32 

5 

2-DEC-1988 

08:32 

4224 

29-NOV-l 988 

11:00 

1075 

30-NOV-l 988 

16:40 

220 

1 L-DEC-1 988 

09:01 

1038 

I2-DEC-1 988 

09:09 

2416 

22-DEC-1988 

08:40 

3 

28-DEC-l 988 

14:25 

2 

1 9- JUL-1988 

22:20 

2 

19-JUL-1988 

22:20 

arnone.lis 
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[DENIS. PCITYJBEND. OUT;  57 
[DENIS . PCITY] BSAS . COM; 6 
[  DENIS .PCITY ] CONTOUR . COM; 3 
[  DENIS .PCITY ] CONTOUR . COM; 2 
[  DENI S . PCITY ] CONTOUR . COM; 1 
[  DENI S . PCITY ] CONTOUR . EXE ; 121 
[DENIS . PCITY] CONTOUR. FOR;  159 
[DENIS . PCITY] CONTOUR. IOP; 9 
[DENIS. PCITY] CONTOUR. LIS; 122 
[  DENIS . PCITY ] CORR . CLN; 2 
[DENIS . PCITY] CORR. CLN; 1 
[  DENI S . PCITY ] CORR . LIS ; 2 
[DENIS. PCITY] CORR. LOG; 4 
[DENIS . PCITY] DAT . DIR; 1 
[DENIS . PCITY . DAT] 19 J1600 . DAT; 1 
[DENIS . PCITY . DAT] 1 9 J1 652 . DAT; 1 
[DENIS . PCITY . DAT] 19 J2022 . DAT; 1 
[DENIS . PCITY . DAT] 19 JUN881 . DAT; 1 
[DENIS . PCITY . DAT] 1 9 JUN882 . DAT; 1 
[DENIS.PCITY.DAT] 20 J1505. DAT; 1 
[DENIS . PCITY . DAT] 20JUN881.DAT; 1 
[DENIS. PCITY. DAT] 21 J1155. DAT; 1 
[DENIS . PCITY . DAT] 21 J1250 . DAT; 1 
[DENIS . PCITY . DAT] 21 J1353 .DAT; 1 
[DENIS. PCITY. DAT] 21 J1441. DAT; 1 
[DENIS.PCITY.DAT] 21 J1519 . DAT; 1 
[DENIS.PCITY.DAT] 21J1 602. DAT;1 
[DENIS . PCITY . DAT] 21 J1630 . DAT; 1 
[DENIS . PCITY . DAT] 21 J1742 . DAT; 1 
[DENIS. PCITY. DAT ] 2 1 J1 8  3 9 . DAT ; 1 
[DENIS . PCITY . DAT] 21 JUN881 . DAT; 1 
[DENIS. PCITY. DAT] 21 JUN882. DAT; 1 
[DENIS.PCITY.DAT] 21JUN883.DAT; 1 
[DENIS . PCITY . DAT] 21 JUN884 . DAT; 1 
[DENIS.PCITY.DAT] 21 JUN885. DAT; 1 
[DENIS . PCITY . DAT] CSRT . DAT; 1 
[DENIS . PCITY . DAT] FOR023 . DAT; 1 
[DENIS . PCITY . DAT] FOR029 . DAT; 1 
[  DENIS . PCITY . DAT ] FRAG . DAT ; 1 
[DENIS . PCITY . DAT] FRAGMENT . DAT; 3 
[DENIS . PCITY . DAT] FRAGMENT . DAT; 2 
[DENIS . PCITY . DAT] IOP020 . DAT;  1 
[DENIS .PCITY . DAT] TRANTEST . DAT; 1 
[DENIS. PCITY ] DECLIST.TXT; 2 
[DENIS . PCITY] DECLI ST . TXT; 1 
[DENIS . PCITY] FOR011 . DAT; 2 
[DENIS. PCITY] FOR0 1 4 . DAT ;  3 
[DENIS . PCITY ] FOR023 . DAT; 27 
[DENIS.PCITYJFOR029.DAT; 3 
[DENIS . PCITY] FOR0 97 . DAT; 65 
[DENIS. PCITY] FRAG. CNV;1 
[  DENIS .PCITY ] GETMCA . COM; 3 
[DENIS. PCITY] IOP020.DAT; 41 
[DENIS. PCITY] LASER. TEK;  20 
[DENIS. PCITY] LASER. TEK;  19 
[DENIS. PCITY] LASER. TEK;  18 
[  DENI S . PCITY ] LN03PLUS . TXT;  1 
[DENIS . PCITY] LPLOT . COM;  3 
[ DENIS . PCITY ] LPLOT . COM;  2 
[DENIS.PCITYlLPLOT.COM; 1 


2 

19-JUL-1988 

22:20 

1 

23-NOV-l 988 

07:42 

1 

20-DEC-l 988 

15:26 

1 

20-DEC-1988 

15:18 

1 

9-DEC-l 988 

08:41 

202 

21-DEC-1988 

14:29 

22 

21-DEC-1988 

11:24 

29 

29-NOV-1988 

08:27 

7  5 

21-DEC-1988 

14:28 

4 

21-DEC-1988 

15:37 

4 

21-DEC-1988 

15:35 

12 

21-DEC-1988 

15:33 

2 

21-DEC-1988 

15:30 

2 

28-NOV-1988 

09:49 

98 

24- JUN-1 988 

14:15 

389 

ll-JUL-1988 

14:40 

170 

24-JUN-1988 

14:36 

586 

25-JUN-1988 

14:03 

272 

24-JUN-1988 

13:55 

134  . 

3-NOV-1988 

09:14 

373 

30-JUN-1988 

11:11 

94 

24- JUN-1988 

14:49 

117 

24-JUN-1988 

14:53 

70 

24-JUN-1988 

14:58 

37 

24-JUN-1988 

15:01 

63 

24-JUN-ly88 

15:02 

36 

24-JUN-1988 

15:05 

101 

24-JUN-1988 

15:06 

38 

24- JUN-1988 

15:11 

88 

24-JUN-1988 

15:12 

169 

24-JUN-1988 

14:34 

116 

24-JUN-1988 

14:59 

206 

24-JUN-1988 

15:15 

211 

18-JUL-1988 

13:01 

175 

18-JUL-1988 

13:16 

1 

3-NOV-1988 

09:13 

2 

4-NOV-1988 

13:35 

1 

4-NOV-1988 

13:50 

6 

7-NOV-1988 

11:57 

139 

3-NOV-l 988 

09:06 

134 

3-NOV-1988 

09:12 

127 

4-NOV-1988 

13:35 

6 

7-NOV-1988 

11:59 

1 

21-DEC-1988 

10:34 

1 

21-DEC-1988 

10:32 

6 

l-DEC-1988 

12:24 

17 

13-DEC-1988 

08:45 

2 

22-DEC-1988 

15:01 

1 

21-DEC-1988 

13:45 

77 

21-DEC-1988 

11:34 

10 

7-NOV-1988 

11:59 

1 

l-DEC-1988 

15:50 

130 

21-DEC-1988 

12:25 

1291 

21 -PEC-1988 

17:26 

1348 

21 -DP P -1988 

17:15 

468 

21-DEP-l 988 

15:10 

1 

20-DEC-l 988 

13:09 

1 

21-DEC-1988 

08:54 

1 

21-DEC-1988 

08:40 

1 

21-DEC-1988 

08:38 

arnone.fis 
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[  DENI S . PCI TY ] Ml . EXE ; 2 
[  DENIS . PCITY] Ml . FOR; 3 
[DENIS . PCITY] Ml . LIS; 2 
[  DENIS . PCITY ] MAK . DAT; 1 
[DENIS . PCITY] NICE1 .OUT; 17 
[DENIS . PCITY ] PCITY . INFO; 4 
[DENIS. PCITY] PCITY. INFO; 3 
[DENIS . PCITY] PCITY. INFO;  2 
[DENIS . PCITY] PCITY . INFO; 1 
[DENIS. PCITY] PCITYMAP.COF;  2 
[DENIS . PCITY] PCITYMAP . DAT;  4 
[  DENIS . PCITY ] PCITYMAP . F10 ; 4 
[ DENIS . PCITY ] PCITYMAP . MCA;  4 
[DENIS . PCITY] PGMS . DIR; 1 
[  DENI S . PCITY . PGMS ] ADDNUM . FOR;  1 
[  DENIS . PCITY . PGMS ] CNVTLT . FOR; 1 
[DENIS . PCITY . PGMS] COLLCT. FOR;  1 
[DENIS . PCITY . PGMS] CONREC . FOR; 1 
[DENIS . PCITY . PGMS] CONTOUR. FOR; 161 
[DENIS . PCITY . PGMS] CONTOUROLD . FOR; 1 
[  DENIS . PCITY . PGMS  ] CONVRT . EXE ; 1 
[DENIS . PCITY . PGMS] CONVRT . FOR; 1 
[ DENI S . PCITY . PGMS ] CONVRT8  8 . EXE ; 1 
[DENIS . PCITY . PGMS] CONVRT88 . FOR; 2 
[DENIS . PCITY . PGMS] CONVRT 8 8 . FOR; 1 
[ DENI S . PCITY . PGMS ] CONVRTEST . FOR; 1 
[DENIS . PCITY . PGMS] CORR. SAS;  1 
[  DENIS . PCITY . PGMS ] GRAPH . FOR;  1 
[  DENI S . PCITY . PGMS ] JULIAN . FOR; 1 
[DENIS. PCITY. PGMS] Ml. FOR;  1 
[ DENI S . PCITY . PGMS ] MERGETT . EXE ; 3 
[ DENI S . PCITY . PGMS ] MERGETT . FOR; 5 
[DENIS. PCITY. PGMS ] MERGETT . FOR;  4 
[DENIS . PCITY . PGMS] NICE1 .OUT; 11 
[  DENI S . PCITY . PGMS ] OPTGRAPH . FOR;  1 
[DENIS . PCITY . PGMS] PARSE . FOR;  1 
[DENIS . PCITY . PGMS] PLOT . SAS;  2 
[  DENI S . PCITY . PGMS ] PLOT . SAS ;  1 
[DENIS . PCITY . PGMS] PRETABLE . EXE; 2 
[DENIS. PCITY. PGMS] PRETABLE. FOR; 3 
[ DENIS . PCITY . PGMS ] PRETABLE . FOR; 2 
[ DENI S . PCITY . PGMS ] PUTALL . SAS ; 1 
[  DENIS . PCITY . PGMS ] READ . SAS ;  1 
[ DENI S . PCITY . PGMS ] READ1 . SAS ; 1 
[DENIS . PCITY . PGMS] READ350 . FOR; 1 
[ DENIS . PCI TY . PGMS ]  READ JUN . EXE ; 1 
[ DENI S . PCITY . PGMS ] READ JUN . FOR; 2 
[  DENI S . PCITY . PGMS ] READ JUN . FOR; 1 
[DENIS . PCITY . PGMS] SET_42 . FOR; 1 
[DENIS. PCITY. PGMS] SUMMARY. SAS; 1 
[ DENIS . PCITY . PGMS] TESSPLOT . FOR; 1 
[DENIS . PCITY . PGMS] TIMETST . EXE;  1 
[DENIS. PCITY. PGMS JTIMETST. FOR;  1 
[DENIS . PCITY . PGMS] TTPLOT . EXE;  5 
[DENIS . PCITY . PGMS] TTPLOT. FOR; 6 
(DENIS. PCITY. PGMS] TTPLOT. FOR;  5 
[  DENI S . PCITY . PGMS ] TUPSPLOT . FOR; 1 
[DENIS . PCITY] PLOT. LIS; 3 
[  DENIS . PCITY ] PLOT . LOG ; 7 
[DENIS . PCITY] PLOTFLUOR. OUT;  1 


43 

2-DEC-1988 

12:30 

12 

2-DEC-l 988 

12:30 

28 

2-DEC-1988 

12:30 

6 

2-DEC-1988 

14:44 

7 

27-DEC-1988 

14:03 

2 

28-NOV-1988 

11:27 

1 

28-NOV-1988 

09:55 

1 

28-NOV-1988 

09:49 

1 

28-NOV-1988 

09:48 

1 

2-DEC-1988 

12:30 

64 

2-DEC-1988 

09:17 

73 

13-DEC-1988 

14:05 

64 

2-DEC-1988 

12:21 

3 

28-NOV-1988 

09:37 

1 

2-DEC-1988 

08:31 

2 

13-AUG-1986 

14:41 

6 

5-FEB-1988 

12:51 

84 

l-DEC-1988 

13:32 

23 

27-DEC-1988 

12:07 

18 

■  8-DEC-1988 

15:19 

20 

28-JUN-1988 

16:01 

21 

5-FEB-1988 

13:01 

22 

12-JUL-1988 

12:09 

37 

27-DEC-1988 

11:52 

37 

12-JUL-1988 

12:08 

22 

2-JUL-1986 

11:28 

2 

27 -DEC-1 988 

13:45 

3 

15-AUG-1986 

16:33 

1 

22-JUN-1985 

13:48 

13 

27-DEC-1988 

13:42 

11 

28-NOV-1988 

11:15 

9 

27-DEC-1988 

11:45 

8 

28-NOV-1988 

11:14 

4 

27-DEC-1988 

14:46 

4 

15-AUG-1986 

16:35 

5 

21-DEC-1984 

08:43 

3 

27-DEC-1988 

14:45 

2 

27-DEC-1988 

13:50 

14 

28-NOV-1988 

11:02 

14 

27-DEC-1988 

12:01 

13 

28-NOV-1988 

11:01 

2 

27-DEC-1988 

13:22 

3 

27-DEC-1988 

13:03 

1 

29-NOV-1988 

15:23 

6 

27-DEC-1988 

14:00 

10 

9-JUL-1988 

11:01 

8 

27-DEC-1988 

11:48 

8 

9-JUL-1988 

11:01 

2 

21-DEC-1988 

13:00 

1 

29-NOV-1988 

14:22 

14 

22-FEB-1988 

16:18 

8 

29-JUN-1988 

08:41 

3 

29-JUN-1988 

09:30 

96 

4-NOV-l 988 

11:43 

13 

27-DEC-1988 

14:40 

12 

22-DEC-1988 

14:05 

19 

2  4 -MAR- 198  8 

16:36 

21 

22-DEC-l 988 

10:51 

4 

22-DEC-1988 

10:50 

53 

4 -NOV-1 988 

11:03 

arnone.lis 
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[DENIS. PCITY] PRETABLE. OUT;  1 
[DENIS . PCITY ] PROGRAMS . COM;  1 
[DENIS.PCITYJPUTALL.LOG;  6 
[DENIS . PCITY ] READ . LOG; 5 
[DENIS . PCITY] READ1 . LOG;  3 
[DENIS . PCITY ] READ 3 50 . FOR; 2 
[DENIS. PCITY] RSAS.COM; 1 
[DENIS. PCITY] SAS.COM; 2 
(DENIS . PCITY] SAS . COM; 1 
[DENIS. PCITY] SET_42 . FOR;  1 
[  DENI S . PCITY ] SET_4  2 . LI S ;  1 
[  DENI S . PCITY ] SUMMARY . LI S ; 1 
[DENIS . PCITY] SUMMARY . LOG;  2 
[DENIS. PCITY]TEMP. LIS;  1 
[DENIS. PCITY]TEMP. MASK;  4 
[  DENI S . PCITY ] TEMP . TEMP ; 2 
[DENIS. PCITY] TEMP. TEMP; 1 
[  DENI S . PCITY ] TTPLOT . EXE ; 2  0 
[  DENI S . PCITY ] TTPLOT . FOR ;  1 9 
[DENIS . PCITY] TTPLOT . INFO; 1 
[  DENIS .PCITY ] TTPLOT . LI S ; 1 6 
[DENIS. PCITY] TTPL0T19. I0P;1 
[  DENI S . PCI TY ] TTPL0T1 9P . IOP ; 2 
[DENIS .PCITY ] TTPLOT20 . IOP; 8 
[  DENIS. PCI TY ] TTPL0T2 OP . IOP; 2 
[DENIS . PCITY] TTPL0T21 . IOP;  1 
[DENIS . PCITY] TTPL0T21P . IOP; 1 


191 

4 -NOV-1 988 

11:04 

1 

25-OCT-l 988 

13:42 

3 

12-DEC-1988 

09:09 

7 

29-NOV-l 988 

10:51 

2 

28-NOV-1988 

13:30 

5 

19-DEC-1988 

13:23 

1 

11 -NOV-1 988 

11 : 49 

1 

30-NOV-1988 

16:06 

0 

30-NOV-1988 

16:05 

2 

21-DEC-1988 

12:59 

7 

21-DEC-l 988 

13:00 

2 

27-DEC-1988 

13:46 

2 

27-DEC-1988 

13:45 

1 

8-DEC-l 988 

13:45 

6 

2-DEC-1988 

17:58 

6 

2-DEC-1988 

17:20 

0 

2-DEC-1988 

17:14 

107 

22-DEC-1988 

14:05 

12 

22-DEC-1988 

13:26 

1 

27-DEC-1988 

14:14 

37 

22-DEC-1988 

14:05 

409 

21-DEC-1988 

17:08 

392 

22-DEC-1988 

14:49 

150 

21-DEC-1988 

15:09 

139 

22-DEC-1988 

14:17 

403 

21-DEC-1988 

17:23 

384 

22-DEC-1988 

15:01 

Total  of  191  files,  33358  blocks 
End  of  save  set 


